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Abstract 

Effective field tlieories provide a formalism for categorizing low-energy effects of a high- 
energy fundamental theory in terms of the low-energy degrees of freedom. This process 
has been well established in mapping the fundamental theory of QCD in terms of the 
hadronic degrees of freedom, which allows for quantitative connections and predictions 
between hardronic observables. 

A more direct approach to performing the non-perturbative QCD calculations is through 
lattice QCD. These computationally intensive calculations approximate continuum physics 
with a discretized lattice to extract hadronic phenomena from first principles. However, as 
in any approximation, there are multiple systematic errors between lattice QCD calculation 
and actual hardronic phenomena. 

To account for these systematic effects in terms of hadronic interactions, effective field 
theory proves to be useful. However, the fundamental theory of interest here is lattice QCD, 
as opposed to the usual continuum QCD. In this work, the basics of this process are outlined, 
and multiple original calculations are presented: effective field theory for anisotropic lattices, 
1=2 vrvr scattering for isotropic, anisotropic, and twisted mass lattices. Additionally, a usage 
of effective field theories and the employment of an isospin chemical potential on the lattice is 
proposed to extract several computationally difficult scattering parameters. Lastly, recently 
proposed local, chiral lattice actions are analyzed in the framework of effective field theory, 
which illuminates various challenges in simulating such actions. 
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Chapter 1 

Introduction 



Quantum chromodynamics (QCD) is the fundamental theory that governs hadrons 
(mesons and baryons) and hadronic interactions. In QCD, fermions in the one-index funda- 
mental representation of the SU{3) gauge group^, also known as quarks, are coupled with 
spin-1 bosons of the same gauge group, known as gluons. This non-abelian SU{3) gauge 
group, whose rank is referred to as color, contains 8 generators, each of which are associated 
with a gluon in the theory. All hadrons are "color neutral"; baryons contain three quarks 
in a state antisymmetric in color, and mesons contain a quark and anti-quark that form a 
color singlet. The QCD Lagrangian including the three lightest quarks (up, down, strange) 
is given by 

1 ^ - 

C = --F-,Fr + Y.^fS]p+m,^)^^^f, (1.1) 

/=i 
where tp is the quark fields with flavor index / and color index i,j, D^^ij = d^Iij + igA^^ij 
with Afj_^ij = A'iTf- and T?- being the eight generators of SU{3), and the kinetic and self 
interactions of the gauge fields are given by F^^ = d^A'j^ — d^A^^^ + ilA^^ , A"]"- . 

The most significant properties of QCD are that the theory is both asymptotically 
free and confining. In other words, at short distance/high-energy, the coupling becomes 
perturbatively small and at long distances/low-energy, the coupling becomes strong. A 
perturbative calculation of the QCD beta function at high energy yields a negative sign (for 
the number of flavors of quarks in the standard model) , which leads to a decaying coupling 
constants as the energy increases (asymptotic freedom). While not mathematically proven, 
QCD is believed to be confining, where it requires more and more energy to separate a quark 
or anti-quark in a hadron. What makes this statement so difficult to prove is the fact that 
at long distances the theory is strongly coupled, and thus, the usual perturbative arguments 
break down. As a result, another method is required to compute non-perturbative physics, 
and lattice QCD is the best option available. 

1.1 Lattice QCD 

In quantum field theory, there exists multiple methods for regulating divergent calcu- 
lations. One such method is to use a lattice; that is, to break up the continuum of space 
and time into many individual sites, which is often referred to as discretizing space and 
time. For every finite sized lattice, there are two parameters introduced: the lattice spacing 



^This is equivalent to two-index antisymmetric representation for an SU{3) gauge group 



between sites and the total length of the lattice. For a simple cubic lattice (hypercubic if 
time is included), these two parameters in a given direction are simply related by the num- 
ber of sites in that direction. By definition, this particular regulator breaks the continuous 
translational symmetry to a discrete translational symmetry, but as we take the continuum 
limit (the spacing between sites goes to zero as the length of the lattice goes to infinity), 
we hope to recover the continuum result. 

Non-perturbative lattice calculations take advantage of this regulator along with the 
path integral representation to perform calculations that do not require a perturbatively 
small parameter. The path integral approach requires defining a partition function given 
by 

Z= [ d[^]d[iP]d[A]e'^'^^'^'^\ (1.2) 



where d[il^] = Yl^ dtpn represents all possible paths of tp and the action is related to the 
Lagrangian by S'(V',^,^) = J d'^xC{ip,ip,A). Another useful quantity to define is the 
correlation function given by 

(/(V^, V,^)> = Z~' /d[V^]d[V']d[A]/(V^, V,^)e*^(^''^'^\ (1.3) 



which is analogous to the definition of the expectation value in statistical mechanics. Both of 
these equations require (infinitely) many multiple-dimensional integrals to compute, which 
is computationally difficult for non-perturbative systems. Rather than the diagrammatic 
approach from perturbation theory, lattice calculations utilize a Monte Carlo integration 
method^, where the exponential of the Euclidean rotated^ path integral acts as a proba- 
bility distribution. As a result of analyzing correlation functions in this framework, energy 
differences and matrix elements can be extracted. The example of the scalar case, as in 
Ref. [2], is shown in Appendix A. 

There are several key issues that arise in non-perturbative lattice calculations. One 
that is relevant to the results in this paper is the fermion doubling problem. As mentioned 
previously, to perform numerical lattice calulations, the space and time directions are dis- 
cretized, with the inverse of the lattice spacing between them acting as the ultraviolet cutoff 
of the theory. However, depending on exactly how the action is discretized (in particular, 
the fermion terms), a single fermion attached to each node can describe multiple fermions 
in the continuum limit. The most basic example of this effect is from the naive fermion 
lattice action. The fermionic part of this discretized action (without gauge interactions) is 
given by 

SnF = «*^ X] ^ [i'x7f^i^x+fi - '^x+illi^i^x , (1-4) 

X,fl 

where, in conventional lattice notation, a is the lattice spacing, x is the location of a given 
site and x+fi represents a site that is one lattice spacing away from x in the /x direction*^, and 
^px and 7^ are the fermionic fields and Dirac matricies, respectively. This action is equivalent 
to taking the continuum kinetic action (tp^f^d^i/j) and approximating the derivative in the 
naive way (dxtpix) ~ ^(^(s^ + a) — 'tp{x)\). Using the discrete Fourier expansion 

p 



^For more detailed information, see Ref. [1]. 
'Transformation defined by going to imagin; 
^H — 1, • • • ,4, where the fourth direction is the Euclidean rotation of the Minkowski time direction 



^Transformation defined by going to imaginary time, t — )• ir. As a result the integrand e — >■ e 



the action becomes 



S^P = T4Y1 V'p7M;r-(e'"^'^ - e-'''P^)i;p = jjY. 'f'p^'^ Iv. Map J V'p, (1-6) 






where ^V'pT^ sin(ap^) is the reciprocal of the propagator. Taking the continuum hmit naively 
(a — 7- 0) yields the correct continuum (Euclidean) behavior for the kinetic term (- sin(ap^) « 
ip^). However, it is important to note that the for every pole of the propagator, there is 
a particle. In other words, a particle is present whenever sin(ap^) = 0. Thus, values of 
ap^ = (0, tt) lead to two poles in each direction. For a four-dimensional lattice, there are 2^ 
poles, which lead to 16 particles that are often referred to as "doublers." In this particular 
discretization, by representing a single fermion field at each node, the action is actually 
describing 16 fermions. Since the goal of lattice QCD is to simulate two or three flavors, 
another discretization is required. 

In order to achieve the correct continuum behavior of the theory being simulated on the 
lattice, an natural extension is to add another term to the naive quark action that vanishes 
in the continuum limit. As a result, the Wilson action was developed, which, in lattice 
terminology is given by 

Sw = a*^ — [V'x(7/. - r)^x+(i - ^x+iiilf^ + r)'>px + rip^^^ , (1.7) 

X,fl 

where r is the Wilson parameter. While this action appears significantly different from the 
naive quark action, it is equivalent to adding the discretized version of the term aripd^d^ip. 
Applying the same Fourier transform as before, 



•S'vK = 74 X] ^p ~^^' sin(ap^) + - ( cos{ap^) - l) 
p,t^ 



V'p- (li 



This action only contains one pole (when ap^ = 0) and reproduces the correct continuum 
behavior as the lattice spacing approaches zero. Thus, this action describes only one particle. 
If we included two or three flavors of this action, we can simulate two or three flavor QCD. 
While this action appears to have removed the fermion doubling problem, it leads to a 
new problem. The issue that arises from this discretization is that it unphysically breaks 
chiral symmetry^. For a free action, this effect will vanish in the continuum limit. How- 
ever, QCD is a gauge theory with gauge interactions that can and will lead to radiative 
corrections. Chiral symmetry restricts many operators that could be radiatively generated. 
However, if chiral symmetry is broken in an unphysical way, as is the case with the Wilson 
action, this will lead to unphysical operators that will contaminate the final results. What 
makes it even worse is that these operators can diverge in the continuum limit and require 
a difficult non-perturbative tuning to remove. These issues will be explained in more detail 
in section 2. 

1.1.1 Chiral Symmetry in Lattice QCD 

One important question in Lattice QCD is whether or not it is possible to both resolve 
the fermion doubling problem and keep chiral symmetry intact. Over the last two decades, 



'The Lagrangian £f = ipOFtp is chirally symmetric if jsDf^s = —Df- A quick way to see if an operator 
is chirally symmetric is to check whether it contains an odd number of 7p's. If it has an even number of 
7p's, it is not chirally symmetric (note, the number of 75 's does NOT alter chirality). 



this topic has been one of great interest and has led to several review articles [3-6] as it is 
now considered a pivotal subject in lattice QCD. 

The Wilson lattice action is an example where chiral symmetry is sacrificed in order 
to remove the doublers. However, this behavior is not unique to the Wilson lattice action 
and, in fact, Nielsen and Ninomiya proposed a no-go theorem for chiral symmetry on the 
lattice. Specifically, the theorem states that a lattice action will contain doublers as long 
as three properties are held true: 1) Discrete translational symmetry, 2) Locality, 3) Chiral 
symmetry. In other words, one or more of these three properties must be sacrificed in order 
to remove doublers. Another way of understanding the Nilsen-Ninomiya theorem is in terms 
of the chiral anomaly. For example, the naive lattice action has an exact U{1)a symmetry 
which persist all the way to the continuum limit. Continuum QCD, on the other hand, 
has an anomolus U{1)a symmetry that is broken by quantum corrections (in particular the 
one loop fermionic triangle diagrams). Thus, in order to produce the correct continuum 
physics for QCD, a source for chiral symmetry breaking must be present to reproduce the 
correct flavor-singlet chiral anomaly defined in the standard way^. One way to do this is to 
break chiral symmetry explicitly, as in the Wilson action. However, as mentioned before, 
this leads to unwanted radial corrections. A more favorable discretization takes advantage 
of the Ginsparg- Wilson relation [8] . 

If a fermionic action has exact chiral symmetry, it obeys the relation 

75D + Z?75 = 0. (1.9) 

Another relation akin to the one above is the Ginsberg- Wilson relation, which is given by 

75Z? + Z?75 = a^75^ , (75^)^ = 75^. (1-10) 

This relation introduces a small source of chiral symmetry breaking at 0{a), and in the 
continuum limit removes doublers, yields radiative corrections ala continuum QCD, and 
yields a source for the flavor-singlet chiral anomaly. Two solutions to the Ginsparg- Wilson 
relation have been found, both of which are equivalent in the appropriate limit. The first 
of such solutions is domain wall fermions [9-11] where four-dimensional chiral fermions 
are given by a mass step function in an additional fifth direction. In the limit of an infi- 
nite fifth direction, this yields a single four-dimensional Dirac fermion (a left-handed and 
right handed "wall" in the fifth direction due to periodic boundary conditions). In reality, 
simulations are performed with a finite fifth directions, but the residual chiral symmetry 
breaking is exponentially suppressed. Another equally valid realization of Ginsparg- Wilson 
fermions are known as overlap fermions [12-14], which satisfy the relations with a non-local 
four-dimensional operator. In has been shown [15,16], that including all the non-local con- 
tributions of the overlap operator is equivalent to the domain wall action with an infinite 
fifth direction (the two actions are related via a KaluzaKlein reduction). 

While both domain wall and overlap lattice actions have favorable features, the major 
drawback to both actions are computational costs. Either simulating an extra dimension 
or multiple non-local operators leads to computational costs about a factor of ten greater 
than four-dimensional, local actions. Thus, it would be beneficial if some cheap, local chiral 
action could be found, and this scenario of minimally doubled, non-orthogonal actions will 
be discussed in chapter 6. 



''As shown in Ref. [7], if one defines a non-singlet axial current for the naive action, the correct chiral 
anomaly can be reproduced. In other words, in the usual notation, the naive action has a flavor non-singlet 
chiral anomaly. 



1.2 Effective Field Theory 

Throughout physics, there are multiple scales that arise when analyzing phenomena. An 
effective field theory (EFT) is a theory that maps the symmetries of a more fundamental 
theory (usually at a higher energy scale) in terms of more relevant degrees of freedom 
(usually at a lower energy scale). An example of this is the low-energy limit of the Standard 
Model, where the heavy, high-energy modes of the W and Z particles are "integrated out" 
and act as small corrections to the low energy theory inversely proportional to powers of 
the W and Z mass. 

To construct an EFT from a more fundamental theory, the first step is to identify 
the symmetries and explicit symmetry breaking in the fundamental theory. Next, for this 
approach to be controlled and fruitful, a seperation of scales is necessary. For example, if a 
theory has a lower scale of p and a higher scale of A, then the expansion parameter in the 
EFT would be p/A. In terms of the relevant degrees of freedom in the EFT, which depend 
on the particular energy ranges of interest, one must write down all the possible operators 
that reflect the symmetries of the fundamental theory (the explicit symmetry breaking can 
be mapped through a spurion analysis) for a given order in the power counting parameter. 
Each of these parameters have a coefficient known as low-energy constants (LECs) that are 
0(c/A") in the power counting where n is determined by the order of the EFT and c is 0(1). 
These LECs are not known a priori and require external input to fix. For non-perturbative 
theories, like QCD, the non-perturbative effects in the calculation are contained in these 
terms, and require a non-perturbative calculation, such as the lattice, to determine. 

Throughout this work, the high-energy, fundamental theory is QCD (more specifically, 
lattice QCD) and the relevant degrees of freedom of interest in our EFT are hadrons, 
nucleons and mesons. The particular EFT of interest will depend which limits we explore. 
For example, if we are looking at characteristic momenta far below the mass of the particles 
of interest, non-relativistic field theory is appropriate. Additionally, if pion mass is small 
compared to the relevant cut-off scales, such as the chiral breaking scale, the pion mass can 
be viewed as a correction to the chiral limit (m^ = 0) . Most of this work will focus on these 
limits, and the resulting EFTs, namely chiral perturbation theory and heavy baryon chiral 
perturbation theory (HB^PT ). 

1.2.1 Continuum Chiral Perturbation Theory 

Chiral perturbation theory (xPT ) is the mesonic EFT of QCD when the characteristic 
momenta are on the same order as the pion mass (771,^ ~ 135 MeV), which are both well 
below the chiral breaking scale (A-^ ~ ^t^"^ f-w ^ 1 GeV). Thus, the characteristic momenta, 
p and the pion mass, m-,^ act as the expansion parameters for continuum xPT . 

As mentioned previously, the first step to constructing xPT is to identify the symmetries 
of QCD. The QCD lagrangian (Eq. (1.1)) has, in addition to the SU{?>) gauge symmetry, 
Lorentz symmetry, and is invariant under parity, charge conjugation, and time reversal. 
However, the main symmetry of focus will be chiral symmetry. In the massless limit, the 
QCD Lagrangian is chirally symmetric, while the mass term breaks the chiral symmetry 
explicitly. This effect can be best seen by rewriting the fermionic part of the QCD lagrangian 
in terms the left and right-handed components of "ip: 

Cp = ^(Jj)^ mq)^ = tpLppL + TpRpl/^R + {ipL'mqi^R + tpRTnglpL) , (1.11) 



where 

ipL = ( — ^)^' ^R=[ — ^)^' 

The chiral symmetric transformation in terms of these components is given by 

where the chiral flavor transformations L and R obey the relation L'L = 1 and R'R = 1. 
The generators of the transformations are those of the SU{Nf) flavor group. In this work, 
we will only focus on either two or three lightest flavors, and these transformations obey 
SU{2) and SU{3) flavor, respectively. 

In the massless limit, the action is invariant under an L and R transformations that are 
independent of each other. As a result, the action pocesses an SU{Np)L0SU{Np)Ri^U{l)B, 
where the last U{1)b represents baryon number conservation. Since QCD is a confining the- 
ory, the chiral condensate ((V'V')) is non-zero, and thus this flavor symmetry is spontaneously 
broken. Both phenomenogical and lattice evidence supports that this symmetry is spon- 
taneously broken to the vector subgroup {SU{Np)l C?> SU{Np)fi — )• SU{Np)v leading to 
Np — 1 Nambu-Goldstone particles. In three flavor QCD, these particles include the pions, 
kaons, and eta (the eta-prime is a flavor singlet whose mass differs from the eta due to 
the chiral anomaly). With the knowledge in hand that non-zero chiral condensate is in the 
vector subgroup {{uu) = {dd) = (ss)), then it is adventagous to write 

ii'LiipRj) = flij = uj5ij, a; / 0. (1.12) 

Performing a chiral transformation, this condensate becomes 

{LnR^)^^=Lo{LR^)^^^Lo^,,. (1.13) 

If L = i?, then Sjj = 6ij leaving the vector subgroup unbroken as expected. If L ^ R, 
then Tiij represents a different vacuum from the vector subgroup. As mentioned before, this 
vacuum, whose generators are that of SU{Np), is a matrix of the Nambu-Golstone particles 
(pion, kaon, eta, etc.) fields, which are the relevant degrees of freedom for xPT . In massless 
QCD, these different vacua are degenerate and each of these particles are massless. However, 
since QCD does contain an explicit chiral breaking term proportional to the quark mass, 
this term "tips" the potential, giving these psuedoscalar particles a mass. To see how this 
occurs, it is important to follow the procedure for mapping the mass term into the chiral 
theory. 

When constructing an effective field theory, one must write down all the possible oper- 
ators at a given order that obey all the symmetries of the fundamental theory. Thus, for 
mesonic observable in xPT , the theory must be written in terms of the mesonic matrix 
Sjj. Since there are no space-time indices in this matrix, the Lorentz symmetry is pre- 
served as long as derivatives of these matrices are contracted in a Lorentz invariant way'^. 



^In Chapter 2 we will explore lattice theories that break the Euclidean equivance of Lorentz symmetry 
and the consequences. As a notational point of emphasis, up and down contracted indices (such as dfj.d'^) 
will represent Minkowski space and only down contracted indices (such as d^d^) will represent Euclidean 
space. 



Under parity, S — )■ S^, and an effective field theory of QCD sliould be invariant under tliis 
transformation^. Under chiral transformations, 



LSi?1" 



RT,^L^. 



(1.14) 



Tlius, wlien constructing tlie cliirally invariant X^T Lagrangian without a mass with the 
monenta p'^ as the expansion parameter, the only combination at leading order (LO) is 
given by 

:^tr(9^St9'^E), (1.15) 

where the coefficient / can be related to the pion decay constant fj,.. Another possible 
operator is tr(S"fS), but since S'I'S = 1, this term has no field dependence is just shifts the 
Lagrangian by an overall constant. Thus, at LO in p^, there is only one term when rUg = 0. 
To include the niq dependence, it is best to rewrite the mass term in QCD in a specific way: 



cTd' = ^tT{d,j:^d^^j:), 



Cm = IpL^qlpR + IpRTult/jL, 



(1.16) 



where mq is the diagonal mass matrix (in the SU{2) isospin limit, niq = ml). This way of 
writing the mass term in QCD seems illogical, but it makes sense when niq is promoted to 
a spurion in order to account for its effects in xPT . Upgrading niq to the spurion^ s + ip, 
this term transforms under chirality the same way as S; namely {s + ip) — )■ L{s + ip)B) . 
Thus, we can now construct the next terms in the chiral Lagangian of the form 



tr((s + ip)T) + (s — ip)^) 



(1.17) 



In the real world and most lattice QCD simulations, niq is purely scalar, so we can set 
s = Bniq and p = 0, which correctly accounts for the chiral symmetry breaking effect of 
niq. B is proportional to the chiral condensate. As a result, all the LO contributions to 
xPT are given by 



40 = :Ltr(a^Sta^S) + ^tr(m,(St + S)). 



The matrix S - 


= ( 


.2*0 is 








<P = 




5 


for two flavors. 
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7r+ 


K^\ 


<P-- 




vr 




2r, 
V6/ 



ni„ 



niu 






rrid 



(1.18) 



(1.19) 



m„ 



TTlu 











md 











m 



(1.20) 



for three flavors. The particular convention used in this Lagrangian is consistent with 
/ ~ 132 MeV. 



**This is not the case for the twisted mass lattice action in Chapter 4 or the convenient lattice action for 
finite isospin density in Chapter 5 

®The s stands for scalar part and p stands for the parity violating pseudoscalar part and should not be 
confused with the momenta. 



As an example of a LO two-flavor tree-level calculation for the pion mass, one must first 
expand (1-18) in terms of pion fields, giving 

4^ = d^TT+d^'TT' - 5o(m„ + md)TT+TT- +■■■ , (1.21) 

where the dots represent three or more pion interactions along with neutral pion interactions. 
From this Lagrangian, the LO pion mass is given by 

ml = Boimu + md), (1.22) 

which is a famous result illustrating that m^ oc niq. Extending this analysis to three-flavor, 
the LO K and r] mass are given by 

m|+ = i?o(m„ + m,) , m^o = Bo{m, + m,) 

2 R /^4 1 1 \ 2 2 ^ 2 2 1 2 n^o, 

m^ = Bo\ -rus + -m„ + -m^ I = ^mj^+ + -m^o - -m^, (1.23) 

where the last equation yields the Gell-Mann-Okubo formula, which agrees with experiment 
up to 10%. The next-to-leading order (NLO) calculation of the mass requires one-loop 
corrections involving the LO terms in the Lagrangian and tree-level contributions from the 
NLO terms in the Lagrangian (terms proportional to m^ or p^). Constructing the NLO 
requires the same procedure as the LO case, namely, writing down all possible independent 
operators that obey the overall symmetry. However, the process is more complicated since 
the number of operators that can be written down has increased signiflcantly. In two-flavor 
xPT , a simplification occurs ((S' + S) happens to be proportional to the identity matrix 
in fiavor space) that leads to only four necessary operators at NLO: 

Ccaat =y tr(a^Sa^St) + :^tr(mgSt + Smg) 

+ j[tr(a^S5^St)]^ + :^tr(a^S9^St)tr(a^St5'^S) 

+ ^^' ^^"^^ [tr(m,St + Sm,)] ' + :^tr(a^Sa^St)tr(m,St + Sm,), (1.24) 

where the LECs ii — ^4 are the original Gasser-Leutwyler coefficients defined in Ref. [17]. 
This Lagrangian will be one of the main focuses in Chapter 3. For three-or more flavors, 
this simplification does not occur and 8 operators are necessary^*^ . The three-flavor NLO 
Lagrangian will not be the focused on in this thesis, but will be included in Appendix B for 
completeness. 

1.2.2 Symanzik Action from Lattice QCD 

As mentioned previously. Lattice QCD is a discretized approximation to QCD. For 
certain lattice actions in the continuum limit (the lattice spacing approaches zero), one 
expects to achieve continuum QCD. However, it is impossible to simulate with zero lattice 
spacings, so it would be convenient to develop a continuum EFT for the lattice action at 
small lattice spacings that would account for both the continuum theory and the flnite a 
corrections that are proportional to powers of the lattice spacings. This EFT is referred to 
as the Symanzik action [18-21]. 



^If the action were to contain external, dyamical gauge fields, there would be more operators needed. 



To illustrate how the Symanzik action is developed, the Wilson lattice action, Eq. (1.7) 
will be used. To begin the process, the symmetries of the Wilson lattice action (now acting 
as our fundamental theory) must be mapped onto the continuum EFT. Many of the sym- 
metries are similar to continuum QCD in Euclidean space, but there are two key differences. 
The first one, that was mentioned previously, is that the action breaks chiral symmetry. 
Second, the lattice does not contain a continuum translation symmetry, but rather contains 
a discretized translational symmetry. To be more specific, the 0(4) rotational symmetry 
of the Euclidean continuum (the Euclidean version of Lorentz symmetry) is broken to the 
hypercubic group. Thus, whatever Symanzik action we write down must account for these 
broken symmetries. 

When constructing an EFT, any possible operator can and will be generated unless a 
symmetry prevents them. As a result, a convenient way of organizing the Symanzik action 
is as follows: 

Hsym = -C^'^ + C^'^ + a£(5) + a2£(6) + . . . , (1.25) 

a 

where a is the lattice spacing, and L'") represents the linear combination of all the linearly 
independent operators with mass dimension n. The dots represents all other powers of 
positive n (there are no operators allowed operators the Wilson action for mass dimensions 
of 1 or 2). Each L^"'' can be written as 

£W=^cJ")0f), (1.26) 

i 
(n) 

with O^ represents the set of allowed independent operators at mass dimension n and 

(n) 

q are the undetermined, unitless coefficients of these operators. As in the continuum 
xPT case, some effort has to be made in order to show that one has the minimum set of 
operators (often, operators can be related through change of variables). Focusing on just 
the fermionic part of the action, the operators needed through 0{a) are given by 






of = ii)a^,uF^,,^lj. (1.27) 



There is only one dimension-5 operator needed, which is referred to as the Sheikholeslami- 
Wohlert [20] term for short. Since there is only one term at 0{a), one approach is to add 
this operator to the actual lattice action to remove all the 0{a) effects in the calculation^^. 
At O(a^), there are 14 additional operators (which are catalogued in Ref. [22]), which 
consist of both bilinears and four-fermion operators. The most notable of these operators 
is ij)^ ^D ^D ^D ^ijj , which is the first term that breaks the 0(4) rotation symmetry to the 
hypercubic group. 

One important issue that becomes clear from the Symanzik action is the complication 
that arise when taking the continuum limit (when a — t- 0). Ultimately, the desired con- 
tinuum limit is QCD. However, as can be seen in Eq. (1.25), there are non-physical terms 
that can either remain or diverge in the continuum limit. There are three categories these 



^^This approach is often referred to as "©(a) improved" or "Clover improved" Wilson action. Since the 
coefficient of the Sheikholeslami-Wohlert term depends on non-perturbative physics, this process requires a 
non-trivial, non-perturabtive fine tuning of the coefflcients in the lattice action. 
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operators fall in: "Relevant" operators that have a mass dimension less than four and di- 
verge in the continuum limit, "Marginal" operators that have a mass dimension of four and 
remain finite in the continuum limit, and "Irrelevant" operators that have a mass dimension 
greater than four and vanish in the continuum limit. For unwanted, unphysical operators, 
irrelevant operators can be suppressed by reducing the lattice spacings. However, relevant 
and marginal operators require an additional non-perturbative fine tuning. For example, 
the Wilson action has one relevant and no marginal operators. The one relevant operator 
~ -ipip has the same form as the QCD mass term. Effectively, this relevant operator acts 
as a large (on the order of the ultraviolet cutoff of the inverse lattice spacings) additive 
renormalization of the quark mass. However, in reality, we want to simulate QCD at finite 
and small quark masses. Thus, when simulating the action, one needs to add an additional 
mass term to the action to cancel this divergence, which requires a non-perturbative tuning 
(such as verifying the measurement of the pion mass in the chiral limit is zero) . A similar 
process would be required for each additional relevant or marginal operator an action might 
generate. It should be noted, however, that each additional fine tuning becomes signifi- 
cantly more difficult to satisfy. Thus, too many fine tunings of a lattice action can prevent 
a lattice action from being practically useful. 

Lattice Chiral Perturbation Theory 

Just as xPT was the EFT of continuum QCD, lattice chiral perturbation theory (L^PT ) 
is the EFT to the continuum Symanzik action in the previous section. The power-counting 
for LxPT that will be analyzed throughout this thesis is p^ ~ vtiq ~ a, known as the generic 
small mass (GSM) power-counting. Another power-counting of interest when performing 
calculations when quark masses are much smaller than lattice spacings, namely ruq ~ o^, 
is known as the Aoki regime. From L^PT , one can calculate the analytic lattice spacing 
dependence of mesonic observables, which aid in continuum extrapolations. This approach 
has been a large focus for categorizing these systematic effects in lattice QCD [22-27]. 

For a Wilson lattice action whose relevant operator has been fine tuned, the new terms 
as compared to continuum QCD are (primarily) the Sheikholeslami-Wohlert term and the 
0(4) breaking term. The 0(4) breaking term, tp^f^D^D^D^ip, first appears at 0{a?'p*), 
which in the CMS power-counting, is at a high order (N^LO) and is a small effect in most 
standard lattice calculations. The greater focus will be on the Sheikholeslami-Wohlert term, 
which obeys the same symmetries as the QCD mass term. Thus, many of the operators 
in LxPT are analogous to the mass terms in xPT with different LECs (indicated by w or 
W). The resulting theory can be given in terms the continuum xPT Lagrangian and the 
contribution from an isotropic lattice spacing. 



AC.SO =^tr(St + S) + ("^+"f-^^° tr(r.,St + Sr.,)tr(St + S) 

+ ^^^^^ [tr(St + S)]^ + ^tr(a,Sa-St)tr(St + S). (1.29) 

This Lagrangian (along with its anisotropic counterpart) and observables calculated is the 
focus of Chapter 3. 
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1.3 Organization 

The structure of this thesis is organized in the following way. Chapter 2, which is 
primarily based on Ref. [28] with P. Bedaque and A. Walker-Loud, develops the EFT for 
anisotropic lattice actions; namely, lattice actions with temporal lattice spacings different 
from spacial lattice spacings. Chapter 3, based on Ref. [29], utilizes this anisotrpic EFT 
to calculate analytic extrapolation formulae for several observables, including pion masses, 
decay constants, and, most notably, phase shifts for 1=2 tttt scattering. Chapter 4, based 
on Ref. [30] with A. Walker-Loud and J-W. Chen, also calculates extrapolation formulae 
for phase shifts for vrvr scattering, but for the twisted mass lattice action, whose results 
were used in practice to extrapolate to the continuum limit in Ref. [31,32]. Chapter 5, 
based on Ref. [33] with P. Bedaque and B. Tiburzi, proposes a novel method of extracting 
several scattering parameters between baryons and mesons, which is currently prohibitively 
expensive in computing time due to the presence of disconnected contributions. Chapter 
6, based on Ref. [34] and previous work in Ref. [35,36] with P. Bedaque, B. Tiburzi, and 
A. Walker-Loud , EFT methods are used to pin down unphysical effects from lattice spacings 
for several non-orthogonal lattice actions, including the recent graphene-inspired lattice 
action proposed by Creutz [37]. As a final note, the reader should be aware that some 
notation has been altered from chapter to chapter to remove ambiguity of symbols (for 
example, in Chapter 4, the lattice spacing is referred to as b as opposed to the usual a since 
a in that chapter is reserved for the scattering length). As a result the conventions are 
defined in each chapter. 
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Chapter 2 

Effective Field Theory for the 
Anisotropic Wilson Lattice Action 

2.1 Overview 

All lattice QCD calculations are performed with finite lattice spacings and finite spacial 
extent (both of which are examples of lattice artifacts). In addition to the extrapolations 
from the unphysically large quark masses that present lattice QCD calculations require, 
extrapolations from a discretized lattice theory to the desired continuum result is also 
necessary. In this work, we will focus mainly on the effects from the finite lattice spacings. 
The formalism that is employed here is a generalization of chiral perturbation theory that 
includes both quark mass and lattice spacing dependence. This patricular formalism, which 
has been utilized to calculate a plethora of lattice observables [22,26,27,38-40], is valid for 
a specific power counting of the lattice spacing effects, specifically these effects are on the 
same order as the contribution from the quark mass (p^ r^ BrUg ~ Wa). These resulting 
formula often depend on LECs which are not known a priori. These LECs, which depend 
on the details of the short distance physics of the non-perturbative calculation, can be 
separated into two categories: ones that contribute to the actual low energy physics of 
interest (physical contributions), and ones that are a result of the finite lattice spacing of 
the particular lattice calculation (unphysical contributions). The goal of this formalism of 
XPT is to isolate these unphysical contributions, so they can be removed in the proper 
continuum extrapolation. 

As mentioned in Appendix A, the lowest energy of the system can be extrapolated from 
the long time behavior of the Euclidean correlation functions. Such correlation functions can 
suffer from signal-to-noise degeneration for certain baryonic observables at long times (see 
Ref. [41] for more details) and can be contaminated by excited states at short times. Thus, 
one method to tame such issues is to utilize an anisotropic lattice where the lattice spacing 
in time differs from the lattice spacing in space by a proportionality factor ^, whose relation 
is given by at = Os/^. Such lattices can give more temportal resolution. This is highly 
beneficial for multiple calculations such as heavy systems consisting of light quarks (for 
example, nucleon-nucleon calculations) where few time slices remain after the excited state 
contamination dies off. In other words, anisotropic lattices can give additional information, 
which might lead to an earlier identification of the ground state plateau. Anisotropic lattices 
have been used extensively in the study of heavy quarks and quarkonia [42,43], glueballs [44, 
45], and excited state baryon spectroscopy [46]. Also, anisotropic lattices may aid in the 
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study of nucleon-nucleon [47, 48] and hyperon-nucleon [49] interactions on the lattice. 

Since anisotropic lattices are in production and currently being used [50-53], the goal 
of this work is to create an effective field theory that encompasses the details and lattice 
artifacts present in these lattices. This work allows for one to perform a systematic analysis 
to identify and remove the new unphysical lattice artifacts that these anisotropic lattices 
present. We begin by constructing the anisotropic version of the Symanzik action [18,19] for 
both the Wilson [54] and 0{a) improved Wilson [20] actions in Sec. 2.2. Then in Sec. 2.3, 
we construct the chiral Lagrangians [17, 55, 56] relevant for these anisotropic lattices for 
both mesons and baryons, focussing on the new effects arising from the anisotropy. We also 
provide extrapolation formulae for these hadrons with their modified dispersion relations. 
In Sec. 2.3.1, we highlight an important feature of anisotropic actions; for a fixed spatial 
lattice spacing and bare fermion mass, if the isotropic action is in the QCD-phase, this does 
not guarantee the anisotropic action is outside the Aoki-phase [57]. 

2.2 Anisotropic Wilson lattice Action 

The starting point for our discussion is the anisotropic lattice action and its symmetries, 
from which we will construct the continuum effective Symanzik action [18,19], which will 
then allow us to construct the low-energy EFT describing the hadronic interactions including 
the dominant lattice spacing artifacts [26]. For the isotropic Wilson (and 0{a) improved [20] 
Wilson) action, this program has been carried out to O(a^) for the mesons [22,26,27] and 
baryons [39,40]. This work is a generalization of the previous work, extending the low-energy 
Wilson EFT to include the dominant lattice artifacts associated with the anisotropy. 

The C'(a)-improved anisotropic lattice action, in terms of dimensionless fields, is given 
by [58, 59] 

S^=Si + Si (2.1) 



t,i<j '^^ t,i n 



atmo + Wt{U) + ^W^{U) 



V'n, (2.2) 



S"! = - Vri 






^n . (2.3) 

. . en 

i<3 

Here, ^o is the bare anisotropy, Pij and Pu are space-space and space-time plaquettes of 
the gauge links V . The bare (dimensionless) quark mass is atmo, Wt{U) and Wi{U) are 
Wilson lattice derivatives and ly is a, parameter which must be tuned to obtain the correct 
the "speed of light." The fields Fti{U) and Fij{U) are lattice equivalents of the gauge field- 
strength tensor in the space-time and space-space directions. S\ is the unimproved Wilson 
action and S2 is the anisotropic generalization of the Sheikholeslami-Wohlert term. The 
coefficients, cj and c^, appearing in S2 are needed for 0{a) improvement of the anisotropic 
Wilson lattice action [20]. At the classical level, they have been determined to be [59] 

C. = F, (2.4) 

where ^ = as/ at is the renormalized anisotropy.^ The choice of using ^ as opposed to ,^0 is 
conventional and the difference amounts to a slightly different value of u. 

^ These parameters can also be tadpole improved with no more effort than in the isotropic action [60]. 
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This anisotropic lattice theory retains all the symmetries of the Wilson action, except 
for the hypercubic invariance, and thus respects parity, time-reversal, translational invari- 
ance, charge-conjugation and cubic invariance. In addition, for suitably tuned bare fermion 
masses, atmo, the theory has an approximate chiral symmetry, SU (Nf) ^i^ SU (Nj) r, which 
spontaneously breaks to the vector subgroup. 

2.2.1 Anisotropic Symanzik Action 

We begin by constructing the Symanzik Lagrangian for the unimproved anisotropic 
Wilson S^ lattice action. This will allow us to set our conventions and introduce a new 
basis of improvement terms which is advantageous to studying the new lattice artifacts 
which are remnants of the anisotropy. In terms of dimensionful fields, the anisotropic 
Symanzik action is given by 

'-Sym — '-Sym ^ "'s^Sym ^ ^s'-Syra ' \'^-'^) 

where we have conventionally chosen to use the spatial lattice spacing as our Symanzik 
expansion parameter. In terms of the dimensionful fermion fields 

1 , 1 



r\j 



the anisotropic Lagrangian is given through 0{a) by 

■^lym = q[lp+mq\q + asq ct auFu + c,. ^ (^ijF'ij Q ■ (2-7) 

i<j 

We have assumed that the parameter v has been tuned in such a way as to make the 
breaking of 0(4) symmetry to vanish in the continuum limit. Otherwise, the quark kinetic 
term would separate into two terms, with an additional free parameter appearing in eq. (2.7). 
To clearly identify the new lattice spacing effects associated with the anisotropy, it is useful 
to work with the basis, 

^Sym = Q W+^q] 1 + "s"^ Q^^t^u^t^u Q + as{ct - Cr) qduFu Q , (2.8) 

from which we recognize the first 0{a) term as the Sheikholeslami-Wohlert [20] term which 
survives the isotropic limit, csw = Cr/2. The second 0{a) term c|^ = (q — c^), is an 
artifact of the anisotropy and the focus of this work. We can classify the effects from this 
operator and subsequent anisotropy operators at higher orders into two categories: those 
which contribute to physical quantities in a fashion similar to the lattice spacing artifacts 
already present, and those which introduce new hypercubic breaking effects. The first type 
of effect will be difficult to distinguish from the already present lattice spacing artifacts 
which survive the isotropic limit. ^ The second category of effects are unique to anisotropic 
actions, and therefore more readily identifiable from correlation functions. 



^Assuming a given set of lattice calculations are close to the continuum limit, and that a range of 
anisotropies is employed, one can disentangle the lattice artifacts associated with the new anisotropic oper- 
ator from those which survive the isotropic limit. 
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A useful manner to quantify these new anisotropic effects is to recognize that the 
anisotropy introduces a direction into the theory, which we can denote with the four-vector 

4 = (1,0). (2.9) 

This ahows us to re-write the anisotropic Lagrangian in a beneficial form for the spurion 
analysis, 

^w = (lW+mq](l + asq csw cr^uF^iy + 4iy^t^w^ o-^a^i.a q ■ (2.10) 

We then promote both asCsw and asCgyyul^ul, to spurion fields, transforming under chiral 
transformations in such a way as to conserve chiral symmetry, 

asCsw — > L{asCsw)R'^ , (asCsw)^ — > R{asCsw)^L^ (2.11) 

O'sclwuj.ui — > L{asSswul,ul)R^ , {asC^gw^i^l)'^ — ^ Ri^sC^sw^i'^i)'^ ^'^ (2-12) 

By constraining agCsw and asCgy^ujjiUu and their hermitian conjugates to be proportional 
the fiavor identity, they both explicitly break the SU{Np)l SU{Np)ji chiral symme- 
try down to the vector subgroup, just as the quark mass term. In addition, we promote 
ttsCg^^u^ui to transform under hypercubic transformations, so as to conserve hypercubic 
symmetry, 

(^scl^ulul — > asc^^^ului A^pAua ■ (2.13) 

By constraining n^ = (1,0), this spurion explicitly breaks the hypercubic symmetry of the 
action down to the cubic sub-group. 

It is worth noting that, in fact, at 0{a), the anisotropic Symanzik action retains an 
accidental 0(3) symmetry in the spatial directions. Close to the continuum, isotropic 
lattice actions retain an accidental Euclidean 0(4) (Lorentz) symmetry as the operators 
required to break this symmetry are of higher dimension and thus become irrelevant in the 
continuum limit [18, 19]. This phenomena is observed in the isotropic limit where the 0(4) 
symmetry is broken by the operator a^qj^Df^Dfj^Dfj^q. For unimproved anisotropic Wilson 
fermions, the breaking of the hypercubic to cubic symmetry (which can be viewed as the 
breaking of the accidental 0(4) to the accidental 0(3) symmetry) occurs one order lower 
in the lattice spacing, at 0{a), and therefore this will likely be a larger lattice artifact than 
the 0(4) breaking of the isotropic action. We now perform a similar analysis for the O(a^) 
Symanzik action. 

0(0^) Symanzik Lagrangian 

In Ref. [22], the complete set of 0{a'^) operators in the isotropic Symanzik action for 
Wilson fermions was enumerated, including the quark bi-linears and four-quark operators. 
From an EFT point of view, it is useful to classify these operators in three categories, 
those operators which do not break any of the continuum (approximate) symmetries, those 
which explicitly break chiral symmetry and those which break Lorentz symmetry. Most 
of the O(a^) operators belong to the first category. Because of their nature, they are the 
most difficult to determine and ultimately lead to a polynomial dependence in the lattice 
spacing of all correlation functions computed on the lattice (which can be parameterized as 
a polynomial dependence in a of all the coefficients of the chiral Lagrangian). The second set 
of operators, those which explicitly break chiral symmetry, can be usefully parameterized 

16 



within an EFT framework as is commonly done with chiral Lagrangians extended to include 
lattice spacing artifacts [22,26,27,39,40]. The last set of operators which break Lorentz 
symmetry, can also be usefully studied in an EFT framework. In the meson Lagrangian, 
these effects are expected to be small as they do not appear until O(p^a^) [22], while in the 
heavy baryon Lagrangian, these effects appear at 0{a?) [40]. To distinguish these Lorentz 
breaking terms from the general lattice spacing artifacts appearing at 0(0^), one must 
study the dispersion relations of the hadrons, and not only their ground states. This is also 
generally true of all the anisotropic lattice artifacts. 

In the construction of the anisotropic action, it is also beneficial to sperate the operators 
into several categories along the lines of those in the isotropic action mentioned above. We 
do not show all of the new operators, as their explicit form will not be needed, but instead 
provide a representative set of the new anisotropic operators which illustrate the new lattice- 
spacing artifacts. In the first category, we begin with operators which in the isotropic limit 
do not break any of the continuum symmetries. Using the notation of Ref. [22], and using 
a superscript-^ to denote the new operators due to the anisotropy, we find for example 

of -^ {of , ^Of } = {qD,]/)D,q, qDtlpDtq] 

of -^ {of , «0f } = {{qi,q){qi,q) , {qitq){qitq)} ■ (2.14) 

In the second category, operators which explicitly break chiral symmetry, we find 

of -^ {of , «0f } = {{qcy^,q){qa^,q) , {qauq){qauq)} , (2.15) 

from which we observe that there is an operator which both breaks chiral and hypercubic 
symmetry. The last category of operators stems from the Lorentz breaking operator in the 
isotropic limit, 

of -^ {of , ^Of} = {q-f^D^D^D^ q , q^iD.DiDi q} , (2.16) 

from which we note that there is only one operator which breaks the accidental 0(3) sym- 
metry down to the cubic group, ^04 , and therefore the dominant 0(3) breaking artifacts, 
in principle, can be completely removed from the theory by studying the dispersion relation 
of only one hadron, for example the pion. Each of these new operators can be written 
in their spuriously hyper cubic-invariant form by making use of the anisotropic vector we 
introduced in Eq. (2.9), 

^Of =uluiqD^IpD,q, 

^Of = u^ul {q-ff,q){q7uq), 

?0f = uf^ui iq(J^xq){qcruxq), 

^Of = ~iAAAa'Q^,'Du'D.>D,, q, (2.17) 



and similarly for the rest of the dimension-6 anisotropic operators, ^0^_g, and ^0\i_ig. In 
this equation, for we have defined 

Most of these operators do not break chiral symmetry, and therefore are present for chirally 
symmetric fermions such as domain- wall [9-11] and overlap [12-14] fermions. We now 
proceed to construct the anisotropic chiral Lagrangian. 
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2.3 Anisotropic Chiral Lagrangian 

Now that we have the complete set of Symanzik operators through ©(a^) relevant to 
the anisotropic Wilson action and the 0{a) improved version thereof, we can construct the 
equivalent operators in the chiral Lagrangian which encode these new anisotropic artifacts. 
We begin with the meson Lagrangian and then move to the heavy baryon Lagrangian. 

2.3.1 Meson Chiral Lagrangian 

We construct the chiral Lagrangian using a spurion analysis of the quark level Lagrangian 
given in Eqs. (2.10) and (2.17). We generally assume a power counting 

niq ~ aA^ , (2.19) 

but work to the leading order necessary to parameterize the dominant artifacts from the 
anisotropy. At LO, the meson Lagrangian is given by^ 

4 = ytr (a^Sa^St) - :^tr (m.BSt + S(m,IB)t^ 

- ^tr ('a.WSt + S(a,W)'f) - ^tr ('a.W^St + S(a,W«)t) . (2.20) 

By taking functional derivatives with respect to the spurion fields in both the quark level 
and chiral level actions, one can show [61] 



\m 



lim^n^, (2.21) 



and similarly the new dimension-full chiral symmetry breaking parameters are defined as. 



' — i_ji yy „r) 



W= hm^ csw '^ %'''^' , (2.22) 



W^ = J™n 4w 



{qcrtiFuq) 



P 

nrr... F..^n^ 

(2.23) 






This anisotropic Lagrangian, Eq. (2.20), has one more operator than in the isotropic 
limit [26,27], which is simply a reflection that there are now two distinct 0{a) operators in 
the Symanzik action. However, for a fixed anisotropy, ^ = as /at, these two 0{a) operators 
are indistinguishable and can be combined into one. Inserting the tree level values of the 
coefficients in the Symanzik action [59], one finds 

W o^csw = v, (2.24) 

W« « c%w = \{^^-^)^ (2.25) 

where the speed-of-light parameter, z/, is determined in the tuning of the anisotropic action. 
The first clear signal of the anisotropy begins at the next order, 0{ap^) for the unimproved 
action and 0{a^p'^) for the improved action. We first discuss the unimproved case. 



We remind the reader we are working in Euclidean spacetime. We are using the convention / ~ 132 MeV. 

18 



These next set of operators are only present for the unimproved action. In the isotropic 
hmit, it was shown there are five additional operators at this order [27] 

2WiiT fa^Ea^Et) tr (a^WE^ + E(asW)t' 



C 



<j>,am 



+ 2H/'5tr ('a^Ea^Et a,WE"f + E(a,W) 

+ 4W6tr ('rngBEt + E(mgB)t') tr ("osWEt + E(a5W)1' 

+ 4W7tr (mgBEt - E(mgB)t) tr (a^WEt - E(a^W)"f 



+ 4W8tr (rngBE+a^WE^ + E(m<^B)tE(asW)n . (2.26) 

For the anisotropic action, there are an additional five operators similar to the above five 
with the replacement of LECs W4_8 — )■ W^_g and a simultaneous replacement of the conden- 
sate W —7- W^. These five new operators, as with the 0{a) operators, are indistinguishable 
from those in Eq. (2.26) at a fixed anisotropy. There are two additional operators at this 
order however, which introduce new effects associated with the anisotropy. 



£ 



W^tT (a^Ea^E^) tr (ului Uw^E^ + E(asW^)"f 



+ W^ti (a^Ea^E"f u^ui kw^Et + s(a,W^) 



(2.27) 



When the anisotropic vectors are set to their constant value, (u^)"^ = (1,0), one sees that 
these operators lead to a modification of the pion (meson) dispersion relation. 



{eI+pD 1 + w 



where 



.asW 

J -K 



.2 , .2^/, , w«^W^^^a,W« 



{Ei+pi)\l + W 



f2 ' f2 



+ ElW^-' 



J TT 



(2.28) 



W = 32(iV/W4 + Ws) , 

W^ = lQ{NfW^ + W|) , (2.29) 

and Nf is the number of fermion flavors. With the 0{a) improved anisotropic action, these 
effects all vanish and the leading lattice artifacts begin at O(a^). The chiral Lagrangian at 
this next order in the isotropic limit was determined in Ref. [22] for which there were three 
new operators. In the anisotropic theory, there are an additional six operators but just as 
with the 0{a) Lagrangian, the three new anisotropic operators can not be distinguished 
from those which survive the isotropic limit unless multiple values of the anisotropy are 
used. The Lagrangian at this order is 






AWk 



a,WtrfE + En - 4IV| a.W^trfE + E^ 
a^W tr (y. - Et) - AW^ a^W^ tr [e - E^^ 
4VF8 (a,W)2 tr (eE + E'^E'f) - 4ll^| (a,W^)^ tr (eE + E+E'f) 



AW!j 



4W| (asW)(asW5) tr f E + E^ 
4VF| (a,W)(asW^) tr (eE + E+E^ 



4H/'|(asW)(asW«) 



tr E 



Et 



(2.30) 
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The last three operators in this Lagrangian vanish for an 0(a)-improved action as they are 
directly proportional to the product of the two 0{a) terms in the Symanzik Lagrangian, 
Eq. (2.10). The other six operators in this Lagrangian receive contributions both from 
products of the 0{a) Symanzik operators as well as from terms in the O(a^) Symanzik 
Lagrangian.^ Therefore they are still present for an 0(a)-improved action but the numerical 
values of their LECs, Wq_^ and Wg_8, will be different in the improved case. 

At 0{a'^p'^), there are nine new operators, three of which survive the isotropic limit, and 
three of which explicitly break the accidental 0(4) symmetry down to 0(3), 



4,av = Qi(«-w)2tr (a^sa^stj + Q2{asW)hT (s^sa^st j tr (s + St 

+ QiiasW^ftv (a^Sa^St) tr (s + St) + Q|(a, W«)2tr (s^Sa^St [s + St 
+ QiiasW^fului tr (a^Sa^St) + Q|(a,W«)\^u^ tr (a^Sa^St) tr (s + S^ 



+ QliasW^fuiui tr (d^^d^^'^ fs + S^ 



(2.31) 



The first operator in each set of three with coefficients Qi, Q\ and Q\ are modifications of 
the LO kinetic operator, while the remaining operators additionally break chiral symmetry. 
The first operator in Eq. (2.31) is an example of the operators mentioned before which 
do not break any of the (approximate) lattice symmetries, and are therefore amount to 
polynomial renormalizations of the continuum LECs. In this example, we see with the 
replacement 

f ^ f + 8Qi{asWf , (2.32) 

the entire effects from this O(a^p^) operator are renormalized away to all orders in the EFT. 
In the anisotropic theory, this also works for the operator with LEG Q\, as this operator 
does not explicitly break 0(4) symmetry. This provides a further example of effects which 
arise because of the anisotropy but contribute to the low-energy dynamics in an isotropic 
fashion. The operators with explicit anisotropic vectors will modify the dispersion relation 
as in Eq. (2.28) but at 0{a'^). For the C'(a)-improved action, these operators will provide 
the dominant modification to the pseudo-Goldstone dispersion relations. 

The last set of operators we wish to address for the meson Lagrangian are those which 
explicitly break the accidental 0(3) symmetry down to the cubic group. There are two 
operators in the meson chiral Lagrangian but they stem from only one quark level operator 
at 0(0^) and therefore all the LO 0(3) breaking effects for all hadrons can be removed with 
the inclusion and tuning of one new operator in the action, ^0\ from Eq. (2.17). The 0(3) 
breaking operators in the meson chiral Lagrangian are 






cHas^^)%A-~^ip~^U tr (s^'sa.stj tr (a.sa^st 
+ cl{asW^f6l/6lJl6f,^ tr (a^,sa,sta,sa.st) , 



(2.33) 



^For conventional reasons [22] we have normalized these operators to the square of the condensates 
appearing at 0{a), but one should not confuse this to mean that these operators vanish for the 0{a)- 
improved action. 
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Figure 2.1: Proposed Aoki phase diagram in Ref. [26] for isotropic lattice action. Region A 
represents continuum-like QCD and region B represents the parity and flavor broken phase. 
Isotropic actions are generally tuned to lie between the first two-fingers of the Aoki-regime 
in the mg — 9^ plane [62]. However, the anisotropic diagram could look quite different than 
this picture, and caution should be taken to ensure the correct phase. 

Aoki Regime 

Aoki first pointed out the possibility that at finite lattice spacing, lattice actions can 
undergo spontaneous symmetry breaking of fiavor and parity in certain regions of parameter 
space [57], as in Fig. 2.1. 

In Ref. [26], Sharpe and Singleton addressed this possibility within an EFT framework by 
extending the meson chiral Lagrangian to include lattice spacing contributions. We begin 
by summarizing the discussion of Sharpe and Singleton which will allow us to highlight 
the new effects that arise from the anisotropy. For clarity of discussion, we consider the 
unimproved two-flavor theory in the isospin limit. Following the notation of Ref. [26], the 
non-kinetic part of the chiral potential can be written 



tr ( S + S1" 



where ci and C2 are functions of the quark mass niq and the lattice spacing a, 

ci~A^(^+gA)+..., 

C2 ~ A^ -| + niqa + a^A^ + . . . . 



(2.34) 



(2.35) 



The 



denote higher order terms in the quark mass and lattice spacing and we are 



ignoring dimensionless numbers of 0(1). Assuming a power counting niq '^ aA^, the con- 
tributions to the vacuum from the C2 term are suppressed compared to the ci term, and 
the vacuum is in the continuum phase with the chiral condensate aligned with unity. The 
Aoki-phase can occur when there is a flne-tuning of the quark mass and lattice spacing 
contributions to ci such that overall size of ci becomes comparable to C2. Parameterizing 
the S-field as 

T. = A + iT-B, (2.36) 
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with the constraint A^ + B^ = 1, the potential, Eq. (2.34) becomes 

V^ = -ciA + C2A^ . (2.37) 

If the minimum of the potential occurs for —1 < ^o < 1? then the vacuum 

So = ( S ) = Ao + ir • Bo , (2.38) 

develops a non-zero value of Bq, spontaneously breaking both parity and the remnant vector- 
chiral symmetry, SU{2)v — > ^(1) [26], giving rise to one massive pseudo-Goldstone pion 
and two massless Goldstone pions. If the minimum of the potential occurs for Amin < — 1 
or Amin > 1 then the vacuum lies along (or opposite) the identitiy, |^o| = 1) |Bo| = with 
the same symmetry breaking pattern as QCD. 

For unimproved Wilson fermions in the isotropic limit, the leading contributions to ci 
are 

ci = f (mgM + asW) . (2.39) 

In the anisotropic theory, there is an additional contribution to the LO potential, such that 

ci ^ cf = ci + fasW^ 



f (nigM + asW + a,W^] . (2.40) 



If the two terms which contribute to W oc 2{qatiFtiq + qCijFijq) are of opposite sign, then 
W^ oc (qcTtiFtiq) may be the dominant lattice spacing contribution to Ci. Therefore, the 
anisotropic theory may be in the Aoki-regime even when the isotropic limit of the theory 
(with the same Ug) is not, and vice versa. This same discussion holds for C>(a)-improved 
actions as well but with a different power counting. 

This analysis carries important consequences for anisotropic actions with domain-wall 
and overlap fermions as well. These actions are generally tuned to lie between the first two- 
fingers of the Aoki-regime in the tuq — g^ plane [62] as in Fig. 2.1. The optimal value of the 
bare fermion mass, which provides the least amount of residual chiral symmetry breaking, 
may be shifted in the anisotropic theory relative to the isotropic value, and furthermore the 
allowed values of the coupling for which there is a QCD-phase may shift as well. 

2.3.2 Heavy Baryon Lagrangian 

We now construct the operators in the baryon Lagrangian which encode the leading lat- 
tice artifacts from the anisotropy. We use the heavy baryon formalism [63,64] and explicitly 
construct the two-fiavor theory in the isospin limit including nucleons, delta-resonances and 
pions. The extension of this to include the octet and decuplet baryons is also possible. This 
construction builds upon previous work in which the heavy baryon Lagrangian has been 
extended in the isotropic limit to include the 0{a) [39] and O(a^) [40] lattice artifacts for 
various baryon observables.^ The Lagrangian is constructed as a perturbative expansion 



■'In those works [39,40], the heavy baryon Lagrangian was extended exphcitly for Wilson fermions and 
in the latter also for a mixed action scheme with chirally symmetric valence fermions. In this work, we 
explicitly address the anisotropic effects for these two scenarios. One can also include the anisotropic effects 
for twisted mass [65] and staggered fermions [66, 67] following the heavy baryon construction in Refs. [68] 
and [69, 70] respectively. However, we are unaware of specific plans to generate anisotropic twisted mass or 
staggered fermions so we do not pursue this here. 
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about the static limit of the baryon, treating it as a heavy, static source as with heavy 
quark effective theory [71,72]. In building the heavy baryon chiral Lagrangian it is useful 
to introduce a new chiral field,® 

a = \/S = exp r^") , (2.41) 

which transforms under chiral rotations as 

a -^ LaV^ = VaH) , (2.42) 

with V defined by Eq. (2.42). One then dresses the baryon fields with a such that under 
chiral transformations, the nucleon and delta fields transform as 

Tijk —^ ViiiVjjiVkk'Ti'j'k' , (2.43) 

where A^ is a two-component flavor field and T is a flavor-symmetric rank-three tensor, 
normalized such that 

Ni =p, N2 = n, 

rin = A++, Tn2 = -^A+, T122 = -^A^, T222 = A". (2.44) 

Both the nucleon and delta are treated as heavy matter fields and are constrained by the 
relations 

N = —^^N, 
2 



2 



T^ = -^T, (2.45) 



where f^ is the four velocity of the baryon which can be chosen in its rest frame to be 
v^ = (1,0). This has the effect of projecting onto the particle component of the field in 
the rest frame of the baryon. The spin-3/2 fields can be described with a Rarita-Schwinger 
field, which in the heavy baryon formalism gives rise to the constraints 

vT = , S-T = 0, (2.46) 

where S^ is a covariant spin vector [63,64]. The simultaneous inclusion of the nucleon 
and delta fields introduces a new parameter into the Lagrangian, the delta-nucleon mass 
splitting in the chiral limit, 

A = mT- rriN , (2.47) 

which is generally counted as A ~ 771,^ in the chiral power counting [63,64,73]. Because this 
mass parameter is a chiral singlet, it leads to a modification of all the LECs in the heavy 
baryon Lagrangian. These particular effects can be systematically accounted for be treating 
all LECs as polynomials in A [74-77]. In the isospin limit, the heavy baryon Lagrangian, 



^This field is generally denoted as ^ but to avoid confusion with the anisotropy parameter, we use a. 
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including the leading quark mass terms, the leading lattice spacing terms and the leading 
baryon-pion couplings is given by 

Cj^T4> = Nv ■ DN + atasNiv ■ u^u^ ■ D)N - 2aMNNtT{M+) 

+ f^[vD + A]T^ + ataET^iv ■ u^u^ ■ D)T^ - 2aMf^T^tT{M+) 

+ 2gANS ■ AN - 2gAAf^S • AT^ + gAN [fl:^'A,f^ej,,Nk + h.c] . 

- 2awNNtr{W+) - 2a«^,iViVtr(W$) - 2awf^T^tr(W+) - 2a«^,i;r^tr(W$) 

(2.48) 

In this Lagrangian, the chiral symmetry breaking spurious are given by 

^A-^- = - ( amj^a + a'mqa' 



wi 



the chiral covariant derivative is 



y (a{W^)^a + a^W^aA , (2.49) 



(Z)^r),,fc = d^'Tijk + v^,%fc + vj;.,r,,vfc + vj:^,T,jk' (2.50) 

and the vector and axial fields are respectively given by 

^M = ^ {^9^^^ - ^^^M^) • (2-51) 

Relative to the Wilson extension of the heavy baryon Lagrangian in the isotropic limit [39,40] 
there is one additional mass operator for the nucleon and delta fields. There are additionally 
extra derivative operators which give rise to the leading modification of the dispersion 
relation for the baryon fields. As in the meson chiral Lagrangian, Eq. (2.20), at fixed 
anisotropy, these two additional mass operators (with coefficients a^^^ and a^^r) can not be 
distinguished from their counterparts which survive the isotropic limit (with coefficients 
aw and aw)- At this order in the lattice spacing, 0{a), the baryon-pion couplings, gA, 
gAN and gAA, are not modified [39]. This Lagrangian, Eq. (2.48) gives rise to the LO 
and NLO mass corrections to the baryon masses. For example, in Fig. 2.2 we display 
the graphs contributing to the nucleon mass at 0{mq), Fig. 2.2(a), 0{a), Fig. 2.2(b) and 

0{mq ), Fig. 2.2(c) and Fig. 2.2(d). The expressions for the quark mass dependence of 
the nucleon and delta masses can be found for example in Ref. [76], and the lattice spacing 
dependent corrections to these masses for Wilson fermions in the isotropic limit can be 
found in Refs. [39,40]. 

At 0{amq), there are two types of mass corrections, those from loop graphs from the 
lattice spacing dependent operators in Eq. (2.48) and tree level terms from operators in 
the 0{amq) Lagrangian. These corrections will be identical in form to those present in the 
isotropic limit, but with different numerical values due to the new anisotropic operators 
which do not explicitly break the hypercubic symmetry. The loop graphs for the nucleon 
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(a) (b) (c) (d) 

Figure 2.2: Diagrams contributing to the nucleon mass at LO ((a) and (b)) and NLO ((c) 
and (d)). Figure (a) is an insertion of the leading quark mass term proportional to aM- 
Figure (b) is an insertion of the lattice spacing terms proportional to aw and o"^. The loop 
graphs arise from the pion- nucleon and pion-nucleon-delta couplings. These loop graphs 
generically scale as m^ but also depend upon A, and away from the continuum limit upon 
the lattice spacing as well. All vertices in these graphs are from the Lagrangian in Eq. (2.48). 

and delta mass corrections are depicted in Ref. [40] in Figs. 1 and 2 respectively. In the 
isospin limit, there are two tree level operators for the nucleon and deltas each, 

^ntI = - nMwNNti{M+MW+) - ni,^NNti{M+Mwi) 

+ tMwf^,T^ti{M+)tT{W+) + ti^^f^T^tv{M+)tv{wi) , (2.52) 

which also function as counterterms for divergences from the above mentioned loop graphs. 
As with the mesons, we must consider the 0(0^) effects to see the first explicit breaking of 
the anisotropy. Even in the isotropic limit, the lattice spacing corrections to the baryons 
at this order proved to be quite interesting, with the presence of the first operators which 
arise from the 0(4) breaking down to the hypercubic group at this order [40]. There are 
several generic O^a^) operators which do not explicitly break the hypercubic group, 

^NTc^ = - nwwNNtr(W+MW+) - nf^^^iViVtr(wi)tr(wi) 
+ twwf^Tf,tT{W+)ti{W+) + 4^^7;r^tr(W$)tr(W$) 
- ni.^NNtT{W+Mwi) + f^^f^T^tiiW+Mwi) , (2.53) 

where the last two operators vanish for the 0{a) improved action and the first two survive 
the isotropic limit. The 0(4) breaking operators appear at this order in the heavy baryon 
Lagrangian because in the heavy baryon theory, there is an additional vector, the four- 
velocity of the baryon one can use to construct operators. For the mesons, the only vector 
is d^ and so the 0(4) breaking operators in the meson Lagrangian don't appear until 
0{a?'p'^) [22]. These heavy baryon 0(4) breaking operators are [40] 

^NT(t> = "s"-4 Nv^v^v^v^N + alti fpVf,Vf,Vf,Vf,Tp . (2.54) 

Similar to these, the 0(3) breaking operators appear at this order. 

However, even though these two sets of operators explicitly break the 0(4) and 0(3) sym- 
metries respectively, they do not lead to modifications of the dispersion relations until higher 
orders [40], and thus function as mass corrections. There is a subtlety of the heavy baryon 
Lagrangian related to reparameterization invariance (RPI) [78] , which constrains coefficients 
of certain operators in the heavy baryon Lagrangian. For example, in the continuum limit, 
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the LO kinetic operator of Eq. (2.48) is related to a higher dimensional operator in such a 
way as to provide the correct dispersion relation,^ 

n2 

C = NvDN — > Nv DN + N—^N, (2.56) 

such that the energy of the non-relativistic nucleon is given by 

For the anisotropic theory, this dispersion relation is then modified by an additional op- 
erator. While in the continuum limit, RPI constrains the coefficient in front of the 1/Mjv 
operator, the anisotropic action gives rise to modifications of this relation. 

- L»2 

C = Nv DN + N — ^A^ — > 
2Mn 

_ /)2 _ 

C^ = Nv DN + N — ^N + ataENv ■ u^u^ ■ DN + as(TKEN{u^ ■ DfN . (2.58) 
2M]\f 

The resulting nucleon dispersion relation from this action yields 



2Mn 



p2 

En{1 + ataE) = Mn{1 + ataE) - 2Nfas{aw + (^w) + WTf- + (^s(^ke {En - M^f ■ 

(2.59) 

For the C'(a)-improved action, the form of this dispersion relation stays the same with 
{atias} —5- {af^al}. Comparing this dispersion relation to that of the pion, Eq. (2.28), 
it is clear that even if the pion dispersion relation were tuned to be continuum like, the 
nucleon dispersion relation would still contain lattice spacing artifacts. This is a simple 
consequence of the fact that the LECs in the heavy baryon Lagrangian (both the physical 
and unphysical) have no relation to those in the meson chiral Lagrangian. 

2.4 Discussion 

In this chapter, the low-energy hadronic effective field theory including the finite lattice 
spacing effects from anisotropic lattice action was developed. This theory allows for one 
to calculate extrapolation formulas for hadronic processes including pions, nucleons, and 
deltas. The main feature that this extension of effective field theory includes is the effects 
of the hypercubic symmetry breaking present in these anisotropic lattices. Additional ex- 
tensions of this effective field theory could be made for mixed action lattices (lattices in 
with different actions for the valence fermions) as in Refs. [79-82]. Additionally, we point 
out that the that these anisotropic lattices can have an altered Aoki-phase as compared to 
the isotropic case. Thus, if the isotropic action is in a QCD phase, the anisotropic theory 
might still be in an Aoki-phase. Overall, this work will aid in the continuum extrapolation 
of correlation functions computed with anisotropic lattices, which will be computed in the 
following chapter for the pion mass, the pion decay constant, and scattering lengths and 
effective ranges for 1=2 vrvr scattering. 



■^In this equation, D\ ^ D"^ - {v ■ Df . 
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Chapter 3 

Isotropic and Anisotropic Lattice 
Spacing Corrections for 1=2 tttt 
Scattering from Effective Field 
Theory 

3.1 Overview 

Numerical scattering calculations in lattice QCD are being performed by several collab- 
orations. These calculations are performed through the analysis of two hadrons in finite 
volume [83-86]. One such scattering that has gained much attention in the field is 1=2 vrvr 
scattering. The two-pion system is both numerically the simplest and the best understood 
theoretically. In fact, in 1966, the scattering of two pions at low energies was calculated 
at leading order in xPT by Weinberg [87]. The subsequent orders in the chiral expansion 
have been worked out next-to-leading order by Gasser and Leutwyler [17] and next-to-next- 
to-leading order [88-90]. Each new order introduces more undetermined LECs. To have 
predictive power, these LECs must be determined either by comparison with experiment or 
lattice QCD calculational results. 

Such numerical lattice QCD calculations (usually involving phase shifts and scattering 
lengths) have been calculated using Wilson lattice actions [47,91-105] along with a several 
other lattice actions [106-109]. Additionally, there are currently only two fully dynamical, 
2+1 flavor calculations of 1=2 tttt scattering, which use mixed lattice actions [110, 111] . 
Regardless of the action, unphysical lattice artifacts due to the finite lattice spacings exist 
in the numerical results of these calculations. Therefore, measures should be taken in order 
to remove these effects so that the results from the lattice can best represent the continuum 
limit. The analysis in this chapter is applicable for both isotropic and anisotropic Wilson 
actions. 

Effective field theory provides a framework by which one can remove these unphysical 
effects. Lattice spacing effects were first made explicit in chiral perturbation theory by 
Sharpe and Singleton [26]. For the Wilson action, the chiral breaking terms that depend on 
the lattice spacing can be accounted for in a similar way to the chiral breaking quark mass. 
Such methods have been extended to mixed-action, partially quenched theories for mesons 
through 0{a) and 0{a'^) [22,27,80], and baryons through 0{a) [39] and ©(a^) [40,82,112]. 
Additionally, Ref. [80] carries out multiple meson scattering calculations (including 1=2 
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TTTT scattering) for mixed lattice actions and shows that for actions with chiral valence 
fermions, mesonic scattering parameters in terms of the lattice-physical parameters will have 
no counter-terms dependent on lattice spacing through next-to-leading order. Alternatively, 
this work calculates these lattice spacing effects^ for 1=2 vrvr scattering for the chiral breaking 
Wilson fermions in both valance and sea sectors, and shows that these finite lattice spacing 
effects first appear in the next-to-leading order counter-terms for this action. 

Many collaborations are now using anisotropic lattices (lattices with different temporal 
and spacial lattice spacings) as opposed the the usual isotropic lattices. Such lattices can 
probe higher energy states (inverse time spacings a^ '^ 6 GeV) and allow for a greater 
resolution (more data points). However, anisotropic lattices lead to new lattice artifacts, 
including terms that explicitly break hypercubic symmetry. Recent work has derived these 
anisotropic lattice artifacts in xPT for 0{a) and 0{a^) for mesons and baryons [28]. There 
are several 1=2 tttt scattering results published for anisotropic Wilson lattices [95,102,105], 
which can benefit from removing these additional lattice artifacts. 

This chapter presents the 1=2 tttt scattering results from the isotropic xPT and the 
anisotropic xPT . The pion mass and decay constant are also determined in this context. 
Sec. 3.2 presents scattering on the lattice and defines the relevant quantity, k cot 6o, used to 
make comparisons between xPT and the actual lattice calculation. Next, in Sec. 3.3.1, the 
continuum scattering theory from xPT is formulated in the context of this work (originally 
worked out before [17,87,90,114]). Then, in Sec. 3.3.2, the continuum result is extended 
for the isotropic Wilson lattice and finally, in Sec. 3.3.3, the result is extended for the 
anisotropic Wilson lattice. 

3.2 Scattering on the Lattice 

The Euclidean two-hadron correlation function in infinite volume gives no information 
about the Minkowski scattering amplitude (except at kinematic thresholds) [115]. However, 
when the correlation functions of two hadrons in a finite box are analyzed, the resulting 
energy levels are given by the sum of the energies of these two hadrons plus an additional 
energy of interaction, AE, which is related to the scattering phase shift, 6i [83-86]. The 
/ subscript here represents the partial wave contribution of the phase shift. In the infinite 
volume, the relation between the total scattering amplitude, T{s,9), and the partial waves 
amplitude, t/(s), is given by 

oo 

T{s,9) = Y,{2l + l)Pi{cos9)ti{s), (3.1) 

1=0 

where s = 4(771,^^ + ^^), and k is the magnitude of the 3- momentum of the incoming particle 
in the center-of-mass frame. The partial scattering amplitude ti{s) is related to the phase 
shift, 6 1 by 



\ s — AmT^"^ 2i \ ^ ~ 4m^^ cot 6i — i 

These relations allow one to compare the calculated scattering amplitude (in xPT ) 
to the lattice calculation of (5;. The s-wave {I = 0) scattering amplitude is the dominant 
contribution to the total scattering amplitude in most low energy scattering processes and 



^ These effects are in addition to the finite volume corrections to 1=2 tttt scattering fi-om Ref. [113] 
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gives the cleanest signal in the lattice calculation. The s-wave projection of the continuum 
scattering amplitude, to{s), is 

1 /-i 
to{s) = - T{s,e)d{cose). (3.3) 

This s-wave scattering amplitude will be the scattering amplitude analyzed throughout 
the rest of this paper. Following the discussion in Ref. [113], through one loop order in 
perturbation theory in Minkowski space, to(s) can be written as 

Jq [S) Zq [S) IIq (s) 

where t^ (s) is the leading order s-wave scattering amplutide, and t^ (s) 

(^0 ' \^)) ^^ ^^^ ^^^^ (imaginary) part of the next-to-leading order s-wave scattering 
amplitude. At this point, it is advantageous to introduce a ii'-matrix, which is defined 
through one loop as 

^(*) - .(LO)f . JNLO,R). s- (^-^^ 

*o («)-*o (s) 

Taking the real part of the reciprocal of Eq. (3.2) and Eq. (3.4) and comparing to 
Eq. (3.5), one gets the relation 



where 



' -Ref ' ]- \r 




Kis) '\tois)J 327rV 
1 \ _ Re to{s) 


COX Oo[Sj^ 


tois)) {Reto{s)f + {lmto{s)f~ 



(3.6) 



It is worth noting that when keeping terms though one loop, Im tQ{s) does not contribute 
(it contributes at the next order). Combining Eq. (3.5), Eq. (3.6), and Eq. (3.7), one arrives 
at the continuum result 



kcot5o(s) = 16TT\fs Re\ — -— ) ss IGtt^/s , , ( 1 

\to{s)J tr Vs)V 



4^^«'^)(s) 



(^) V i'^'^Hs) 



{3.i 



As previously mentioned, lattice scattering calculations are performed in Euclidean space 
at finite volume. The Euclidean amputated four-point correlator from the lattice, tq{s), is 
given by 

"° ' " tr'(») - tr°'''\') - A<o(») - ig^5(fcig5ni:) ■ 

where Ato{s) represents all of the non-physical lattice artifacts (lattice spacing errors, finite 
volume errors ,etc.), s is related to the energy of interaction, AE, and 5 is a universal 
function of s [85, 116, 117]. If both pious in the box start with no external momentum, then 
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1 

K{s) 


Ato(s) 


1 c(is 
167r\/i 


167r2 







k cot (5( 




1 qf{s- 


-4m,r" 


^) 



s = (AE' + 2m^)^. In this paper, the only effect from lattice artifacts that will be included 
in Ato is the lattice spacing effect. Manipulating Eq. (3.9): 



(3.10) 



The energy states are given by the poles of Eq. (3.10), which are given by [113] 

k cot 5o + A(fc cot 6o) = :^s( ^' ~^^2 ^^ ) ' (3-11) 

where 

A{k cot 6o) = -lQ7^^/^—^^. (3.12) 

In general, most lattice calculations give their results in terms of the scattering length, 

a^^^. One can extract the scattering length and the effective range, r^^^, via the expansion 

of A; cot 6o: 

1 1 

— 5 

2 

It is important to note that the prescription given above for finding the scattering 
length and effective range implies that the lattice artifacts, A(fccot (5o), have already been 
subtracted before the expansion. In this paper, continuum results are given in terms of 
kcotSo and lattice artifacts are given in terms of A(kcot6o). 

For results given in terms of the scattering length and effective range, one can relate 
Eq. (3.13) to the left hand side of Eq. (3.11) to arrive at 



kcot6o = ^ + ^rifk^ + • • • . (3.13) 



fecot^o + A(fecot(5o) = ( i + a(^)^ + i( rif + Ar^f ) fc^ + • • • , (3.14) 



where 



and 



1=2 ' ^\I=2J I ■ 2 



A -^) =A(A:cot5o)|fe2=o, (3.15) 



1=2 _ ^d{A{k cot 5o)) 



^'^^'^ "^ dk^ 



(3.16) 

fc2=0 

While these relations are not too complicated, they do add additional steps to the 
calculation when compared to working with only k cot Sq and A(/c cot 5q). Therefore, if one 
wants to extract the scattering length and effective range from the lattice calculation, one 
should first subtract A(/ccot5o) from the right hand side of Eq. (3.11) and then expand 
to determine the individual parameters^. This paper relates kcot5o and A(kcot6o) to the 
effective field theory of the lattice. 



^Current numerical calculations can only determine fccot(5o for a limited number of k values. This leads 
to inaccuracies in expansions of k^ and adds difficulty to finding the effective range. 

30 



3.3 Chiral Perturbation Theory Results for /ccot^o and 

A{kcot6o) 

In leading order and next-to-leading-order chiral perturbation theory, it is necessary to 
introduce several undetermined LECs in order properly account for corrections and counter- 
terms. The number of independent LECs in the continuum depends on whether there 
are two flavors or more. For 1=2 tttt scattering being calculated here, only two flavor 
xPT {SU{2)l SU{2)ji chiral symmetry) is needed for extrapolation. 

3.3.1 Continuum 

From the continuum xPT Lagrangian, one can predict numerous results for different low 
energy processes involving hadrons. However, the extent of the accuracy of these predictions 
are ultimately tied to how well the LECs are known. For this reason, there has been much 
effort in the lattice community to try to determine these values [118-122]. 

The continuum Lagrangian in xPT is determined order by order in Bruq and p'^. The 
Lagrangian through 0{p'^) for two flavors is given by [17] 



Ccont =4rtr(5^Sa^St) + ^tr(mgSt + Sr 



8 



+ -^[tr(a^sa^st)]^ + :^tr(a^s9^st)tr(a^sta^s) 



4 



+ 



^^' ^ .^^^^ [tr(m,St + Sm,)] ' + ^trid^W^^Mmg^^ + Sm,), (3.17) 



where / ~ 132 MeV, ^i_4 are the original Gasser-Leutwyler coefficients defined in Ref. [17] 
and 



^-M?f), *-(;^:^j. -",-(::)■ (3.i«) 



At LO, the resulting condensate is 

B = lim K|^. (3.19) 

From Eq. (3.17), one can calculate the physical values for the mass of the pion (to-tt) 
and the pion decay constant (/tt) to LO and NLO. These expressions are given by 

m^2 = rr? ^—^\\m^^iX{m^)-m^iX{m^)\^Uz^ (3.20) 



-^il{m^) + 2^4 — 



fn = f l-^r^il{m^) + 2i^-^ (3.21) 



where m? and X{mT^) are defined below in Eq. (3.23). When evaluating the scattering 
amplitude from xPT , one has the option of either expressing the answer in terms of 
the bare parameters (/ and m) or in terms of lattice-physical parameters (/vr and m.,^). 
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(a) 




(d) 





Figure 3.1: The one- loop diagrams which contributing to the rnr scattering amphtude (same 
diagrams that contribute to the finite volume effects in Ref. [113]). The top three represent 
the s-channel (a), t-channel (b), and u-channel diagrams (c). The diagram resulting from 
the six-pion vertex is given in (d) and the wavefunction corrections are in diagram (e). 
These are the only diagrams for continuum, isotropic, and anisotrpoic Wilson actions. 

Throughout this paper, the bare parameters will always be eliminated from the scatting 
amplitude. 

The continuum 1=2 vrvr scatting length at infinite volume is given by 

2 r 2 2(3s-4m^2) {s-2m^'^f ^, 
Tcont = — T^\s- ^rn^ --p2 tl{m^) H -^ ij{m^,ps) 
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'ti (27r)4 fc2 _ ^2 

(fk 1 



In (27r)4 [{k + Pf - m2] [P - m2] ' ^^'^^^ 

This scatting amplitude includes all the partial wave contributions (this T is the same as 
the T{s,6) that appears in Eq. (3.1)). When projecting on the s-wave, expanding through 
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©(fc^/rriTr^) and using Eq. (3.8), the result for /ccot(5o(for the regularization and renormal- 
ization scheme defined in Ref. [90]) is 



k cot 60 



8ttU' 



m. 



m^ 



(47r/^) 



31n(^)-l + ^JJ^) 



2 V (47r/^)2 L 3 



17, /m^^\ 31 ^„ , , 



+ 



m^ 



(3.24) 



where ^"^(/^) and ^^^^(/u) are Hnear combinations of the Gasser-Leutwyler coefficients 
given by [90] 



C(m) = - ^47^f{4{£f{^,)+£f{^,)) + (^f (^) -^f (^))), 



C(m) =4(47r)2(l2^f (^) + 4^f (^) + 7^f (^) - 3^f (^)). 



(3.25) 



The superscript R represents the renormahzed Gasser-Leutwyler coefficients with scale- 
dependence. 

In order for these continuum predictions from xPT to be useful in a physical context, 
one needs to determine the values for i'^.,^{fi) and ^^^{fj,), which are undetermined from 
xPT alone. While numerous values have been quoted for the Gasser-Leutwyler coefficients 
[122], it is still beneficial to determine these values with more precesion. Therefore, it is 
prudent to use various lattice calculations at different pion masses determine these values. 
However, since the lattice calculations are performed with discretized space and time, one 
needs to remove these lattice artifacts to extract the continuum result. 

3.3.2 Isotropic Discretization 

To calculate finite lattice spacing corrections to 1=2 vrvr scattering for the isotropic 
Wilson action, one can follow the same steps done in the continuum case, but starting from 
a Lagrangian which includes these lattice spacing artifacts. The analysis on lattice spacing 
effects was done for the Symanzik action by Sheikholeslami and Wohlert [20]. From this 
analysis, the Lagrangian was made explicit in xPT by Sharpe and Singleton [26], followed 
by Bar, Rupak, and Shoresh [22,27]. The power counting they used for this xPT Lagrangian 
is 



a,W 



Bniq ^ 



P^ ~ e, 



(3.26) 



where a^ is the lattice spacing (same spacing in space and time direction) and VF is a 
condensate defined below. The simplified two fiavor Lagrangian to 0{e'^) is 



'~'iso — '~'Cont ~r ^l-iso' 



(3.27) 



A£.. =^tr(St + S) + ("^ + "f-^^° tr(m,St + Sr.,)tr(St + S) 

+ !%E)! [tr(St + S)] ^ + ^^tr(5,Sa'^St)tr(St + S) (3.28) 
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This Lagrangian is similar to Eq. (3.17) with one new term at LO and three new terms 
at NLO. An new terms are proportional to a^ or a^, which will vanish in the continuum 
limit as a^ — )■ 0. At LO, there is a new condensate given by 

W= hm csw ^ • (3.29) 

The new LECs at NLO are given by w^, tf3, and W4. All these new terms obey the same 
symmetries as the Lagrangian in the continuum case and break chiral symmetry in a similar 
way to the quark mass. It should also be noted that these new LECs depend on a^ In Og as 
well (as opposed to the Gasser-Leutwyler coefficients that have no dependence on the mass 
term) . 

Furthermore, Ref. [123] showed that the axial current (needed to calculate f-,^) has an 
additional term at this order in xPT given by 

AA;^ = 2awAdf,iv{a^{T. - St)) . (3.30) 

This term (which can also be inferred from Ref. [124]) leads to modifications of the LECs 
as well as the coefficient in front of the chiral logarithm in /t^. Thus, this coefficient has 
dependence on the lattice artifacts at NNLO. Ref. [123] also points out that the condition 
for fixing the renormalization factor, Za, of the lattice currents needs to be mapped onto 
xPT ^. From this Lagrangian, /^ and m^ through NLO are [22,27,123] 

mj^ = {w? + 2asW) + —-r[Am^'^iZ{'m^) - (m^ + 2asW)iX{'m^)] 

3/^ 

+44y^ + 8u;3 — -^— + I6W3 ^ (3.31) 

r 2 m'^ a Wi 

U = f l--^il{m^) + 2i^-^ + Aw,ff^ , (3.32) 

where tfe// in fir includes W4^ and wa and can vary based on the given renormalization 
condition for the axial current. 

To acquire the continuum result of these quantities, one needs to remove all the terms 
with dependence on a^ or a^. The resulting 1=2 vrvr scatting amplitude with the physical 
parameters restored is given by 

J-iso ^ J-cont + ^-'-iso-i yo.oo) 

where the m-,^ and f-,^ in Tcont are given by Eq. (3.31) and Eq. (3.32) and ATiso is given by 



AT,3o = -^\- W{W3 - 2^3) "'^7" - 32(w'3 -m + h) 



2 



W 



^2/ asWXs-2m^. 

+ S{weff-U) Tl. \- (3-34) 

Jit 



'^Thc condition for fixing Za is chosen by individual lattice calculations. Ref. [123] shows an example 
of this in xPT by using the chiral Ward identities in infinite volume, which leads to /,r being free lattice 
artifacts until NNLO. 
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It is worth noting that by restoring the physical parameters, m^ does not appear, and 
ttgW only appears in the terms containing the LECs (the i and w terms). This is a bit 
different than the continuum case where one could eliminate m^ with only ttt-tt^- Now, one 
eliminates m? with (m^^ — 2asW), and thus, several LECs are multiplied by factors of UsW. 
In addition, all of the continuum results without LECs remain unchanged since each vertex 
will only contribute 171,7^^ when the physical parameters are restored. Using the relation in 
Eq. (3.12), the resulting artifact for the isotropic Wilson lattice, A(/ccot (5o)iso is given by 






A(fccot5o)..o « ^ <J^)-^ +<^(m) 



-lwl^{fi)^ + 7w':^{^y-^^)^ + --- \, 3.35 



where 



<.(/") = 


- Siinfiw^ifi) - w^ffif,) - 2i§i^^)+ifi^^)), 


<.(m) = 


-n47rf{wf{f,)-w^{f,)+i§{fi)), 


<.(m) = 


- 8(4^)2 (7u;3«(^) - Sw^ffifi) - Ui§{i,) + 3£f (^)) 



(3.36) 

As seen in the results, the artifacts from the isotropic lattice that are present in the 
final form are either linear or quadratic in Ug- By using results that differ in lattice spacing, 
one can pick off the coefficients of these artifacts and remove them from the final result. If 
one is working with a perfectly clover-improved Wilson lattice, this would remove all 0{as) 
effects leaving only the O^a^) effects. It is also important to note that there is no physical 
information gained by determining specific LECs that are a result of the isotropic lattice 
spacing (the individual w terms) unlike determining specific Gasser-Leutwyler coefficents. 
Therefore, the useful coefficient to extract is the linear combination of these terms so they 
can be removed from the final result. 

At this point, one can compare these lattice spacing effects for the Wilson action to 
those found for the mixed action (with chiral valence fermions) in Ref. [79,80]. For this 
mixed action case, when in terms of the lattice-physical parameters, there is no lattice 
spacing dependence through the NLO counter-terms. In contrast, when the lattice-physical 
parameters are restored in the Wilson action, these effects first appear at the NLO counter- 
terms. Thus, these additional effects that are not present at NLO in the mixed action 
calculation will need to be removed for the Wilson action calculation. 

3.3.3 Anisotropic Discretization 

With the recent formulation of xPT for the anisotropic lattice [28], one can begin 
calculating corrections to various quantities of interest on the lattice. This process is, in 
general, very similar to calculations in the continuum and isotropic lattices, but one now 
picks up additional terms that are a result of having different spacial and temporal spacings. 
To help extract these effects in a more simplistic notation, the superscript ^ has been added 
to all the new terms resulting from this anisotropy. In practice, the anisotropic lattice picks 
up two new non-perturbative parameters: the parameter S, = as /at which is the measure of 
anisotropy and the parameter v, which is used to correct the "speed of light" [58,59,125]. 
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By setting both parameters to 1, the isotropic hmit is recovered. In addition to the W 
condensate defined in Eq. (3.29), we pick up a W^ condensate that is given by 

W^= Hm 4^(n«)^n«).i^^^^^. (3.37) 

where n^ is a vector that breaks hypercubic invariance. It is important to note that the con- 
vention chosen for u^ appears in the anisotropic xPT and its observables. For convenience, 
we choose this vector to be n^ = (1,0). The condensates and the anisotropic paramaters 
are related at the classic level by (with Wilson coefficients r^ = r^ = 1) 

W oc csw oc v, (3.38) 



W^ c 4^^-{-l-u]. (3.39) 



2 Us 

In the isotropic limit when ^ and v are set to 1, the isotropic condensate will remain and 
the anisotropic condensate will vanish. Using a similar notation throughout, all terms that 
appear with a ^ superscript will vanish when Ug = at and z/ = 1. 

The power counting convention used in Eq. (3.26) for the anisotropic Lagrangian is 

QsW ~ asW^ r^ Bmq r^ p^ r^ e. (3.40) 

Writing this new Lagrangian in the form of Eq. (3.27), the two- flavor anisotropic 
XPT Lagrangian through ©(e^) is 

'~'aniso — J~-cont ~T ^'-"iso ~r ^J~-aniso- V'^'^-'-j 

where 

A£.„.„ =^^tr(St + S) + M±^i±^4KE!A)^^(^^^t + Sr.,)tr(St + S) 
+ ^'^^f^' [tr(Et + E)]- + 4iasWXasW^) ^^^^^, ^ ^^^^ 

+ !^;4^^tr(a^Sa^St)tr(St + S) + !^il^^^A*u,(5^S5'^St)tr(St + S). 

(3.42) 

In addition to the anisotropic condensate W^ mentioned above at LO, there are five 
new LECs at NLO as a result of this anisotropy. Four of these new LECs obey the same 
symmetry structure as the isotropic terms, however the w^ term additionally breaks hyper- 
cubic invariance. Therefore, this term only corrects the time derivative, but not the spacial 
one (for the convention of uji chosen). As a result, /^^ is parameterized by two constants; 
/*, which is /tt measured in time, and /*, which is f-j^ measured in space. This leads to one 
correction for the space- measured /^ and a separate correction for the /^. The pion mass 
(m-Tr), the time-measured pion decay constant (/^), and the space-measured pion decay 
constant (/^) through NLO are 
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(3.43) 
(3.44) 
(3.45) 



where, as in the isotropic case, the w^^r depends on the renormahzation condition for the 
axial current (this LEG is the same for both /^ and /*). 

For the rest of this section, all calculations are given in terms of m^ and /*. As mentioned 
before, how one accounts for the effect of this hypercubic breaking term depends on the 
convention. In the convention used here, only /^ sees the effect of this term and /* does 
not. The scatting amplitude is given by 



'-cont 



+ ATiso + AT„, 



(3.46) 



where the m^ and /^ in Tcont ai'e given by Eq. (3.43) and Eq. (3.44) and ATaniso is given 

by 



AT,, 



in 
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lQ{wl - 2^3) ^^^)2 32(4 -wl+ £3) ^^,^2 



32{wl -W3-WI + 21. 



{asW){asW^) 



+ ^{wlff + w\-£i) 



ifi? 

asW^js - 2m^^) 



16w' 



t\2 



' (/*) 



(3.47) 



Most of the effects in this scattering amplitude are similar to the isotropic case, except 
now there are also expansions in agW^ in addition to the expansions in agW. Thus, as 
expected, if all the anisotropic effects are removed, only the isotropic limit remains. The 
only new symmetry breaking effect is the wi term which is not a hypercubic invariant 
term. In other words, all of the hypercubic breaking due to anisotropy at this order is 
contained in this term. However, it's effects in ATaniso appear as just another contribution 
to the linear combination of the LECs in front of the term asW^. Therefore, it is difficult 
to determine the effect of the the hypercubic breaking term alone from ATaniso since its 
effects will be mixed in with the other anisotropic LECs*^. The resulting artifacts for the 
anisotropic Wilson lattice are the isotropic artifacts, A(fccot 5o)jso) and the anisotropic 
artifacts, A(A;cot (5o)araso- Therefore, the total effect of the lattice artifacts due to lattice 
spacings are 



The total hypercubic breaking effect would be more visible from the differences of /^ and /^ 
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A{k cot 5o) = A{k cot (5o)jso + A(/c cot 5o)aniso- (3.48) 

The anisotropic lattice artifacts are given by 

+ •••1 (3.49) 

where 



t/;rj/^) = - 8(47r)2(u;f (m) - w%{^i) + t.f (m) - 2£3«(m) + ^f (m)), 



W 



TTTT 



(^) = - 16(4^)2 {wfiii) - w^iii) + i§{f^)) , 



w 



(^) = - 16(47r)2(u)r(^) - wiif,) - t.r(^) + 2i§{n)), 



wilifi) = - 8(47r)2(7«;f (^) + dw^ffif,) + u;f (^) - Ue§{f,) - <(^)). (3.50) 

As a result, the ultimate effects of the anisotropic lattice on A{k cot 5q) are more terms 
that will require variation of Ug and at independently to fit. Three more terms require fitting 
to correct the constant term in the expansion and one more term requires fitting for the 
©(/c^/mTT^) term in order to remove it's effects. As mentioned with the isotropic correction, 
no physical information is gained by picking off the individual anisotropic LECs. Analogous 
to the isotropic case, the anisotropic Wilson action first has lattice spacing dependence at 
the NLO counter-terms. The aim here is to determine these linear combinations and remove 
their effects from the lattice measurements. 

3.4 Discussion 

In this work, 1=2 tttt scattering was calculated from the xPT for isotropic and anisotropic 
lattice spacings. Also, connections between these xPT calculations and the kcotdo value 
measured from lattice calculations were illustrated. When A(/c cot 5o) is given in terms of the 
lattice-physical parameters, these lattice spacing effects first appear at the NLO LECs and 
can be removed from the result of the lattice calculation. However, A(A; cot 6o) has numerous 
undetermined linear combinations of LECs, which would need to be determined (by fitting 
several different lattice spacings) in order to successfully remove it from the lattice result. 
Therefore, as more lattice calculations of tttt scattering are completed (for both isotropic 
and anisotropic lattice spacings), these combinations of LECs can be determined better, 
which will result in a more accurate result after these lattice artifacts are removed. 
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Chapter 4 

TTTT Scattering in Twisted Mass 
Chiral Perturbation Theory 

4.1 Overview 

Calculations of hadron interactions from lattice QCD have received significant attention 
in recent years [126], with dynamical calculations of two-meson systems [101,110,111,127, 
128], two-baryon systems [48,49] and systems of up to 12 pions [129,130] and kaons [131]. 
Additionally, lattice methods are currently being applied to the low-energy effective field 
theory of multinucleon interactions [132-139], for which there exists a nice review [140]. 
The majority of the dynamical lattice QCD calculations of hadron interactions to date, 
have either been performed with Wilson fermions or a mixed lattice action [141, 142] of 
domain-wall valence fermions [9-11] and the Asqtad improved [143, 144] rooted staggered 
MILC configurations [145,146]. Twisted mass lattice QCD [65,147] has recently emerged 
as a viable fermion discretization method for lattice calculations with two flavors of light 
quarks, {up and down) [119, 148, 149] and hopeful prospects of 2 -|- 1 -|- 1 {up, down, 
strange and charm) flavors of dynamical sea fermions [150] in the chiral regime. Therefore, 
it is natural that hadron interactions will be computed with the twisted mass fermion 
discretization method, and initial twisted mass results for 1=2 vrvr interactions have been 
recently published [31,32]^. 

To truly understand the twisted mass results of 1=2 Trvr interactions, the systematic 
effects due to lattice artifacts must first be understood. Lattice discretization effects can be 
incorporated into the chiral Lagrangian through a two-step process first detailed in Ref. [26]. 
One first constructs the effective continuum Symanzik Lagrangian [18, 19] for a given lattice 
action. One then builds the low energy chiral Lagrangian from the Symanzik theory, giving 
rise to new unphysical operators with their own LECs. These new operators capture the 
discretization effects for a given lattice action. This process is spelled out for the Wilson 
action in Chapter 1. 

In this chapter, we briefiy review the construction of the twisted mass chiral Lagrangian 
in Sec. 4.2. We then determine the lattice spacing (6) corrections to low-energy tttt scattering 
specific to the twisted mass lattice action. We work through C(6^) ~ 0{bm'^). 



^ These references made use of the work developed in this chapter 
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4.2 Twisted mass lattice QCD and the continuum effective 
action 

The twisted mass chiral Lagrangian was determined previously in Refs. [124,151-154], 
and for baryons in Ref. [68]. In this work, we focus on twisted mass lattice QCD with 
degenerate light flavors given by the lattice action 



^ = E ^(^) \l E ^'^(^^ + ^'^) - ^ E ^^^- + 



rriQ + i-i^nno 



^(x), (4.1) 



where ij^ and V' are the dimensionless lattice fermion fields, Vi/(V^) are the covariant forward 
(backward) lattice derivatives in the v direction, mo is the dimensionless bare quark mass 
and ^0 is the dimensionless bare twisted quark mass. The fermion fields are flavor doublets, 
Ts is the third Pauli-spin matrix and the bare mass term is implicitly accompanied by a flavor 
identity matrix. Our twisted mass xPT analysis also holds for dynamical lattice calculations 
with 2 + 1 + 1 flavors, the only difference being the numerical values of the LECs determined 
when fitting the extrapolation formula to the calculation results. This action is beneficial in 
several ways. While it has a term that breaks parity, it contains a parity (^ flavor symmetry 
that prevents unwanted radiative corrections. However, this term, which is referred to as the 
twisted mass term, is a mass term that multiplicatively renormalized akin to how the quark 
mass term is in continuum QCD. This prevents issues involving exceptional configurations. 
Additionally, for certain choices of rriQ and /Uq parameters, known as maximal twist, this 
action can additionally be 0{a) improved without including a Sheikholeslami-Wohlert term. 
The continuum chiral Lagrangian, supplemented by discretization effects is determined 
with the two step procedure of Ref. [26] . This was done for the twisted mass lattice action 
in Ref. [153], to NLO in which a power counting niq ^ bAQQj^ was used and which we shall 
adopt. The resulting effective Lagrangian is 

Ceff = ^giue + q{P+ m + 175X3^)^' + csw^ qicTiMuF^uq + 0{b'^, brug, m^) , (4.2) 

where Cgiue is the Yang-Mills Lagrangian. The quark fields are an isodoublet, q^ = {qu, Qd) 
and the quark masses are given by 

m = Zmimo - mc)/b, 
H = Z^Ho/b. (4.3) 

The symmetry properties of the twisted mass lattice action protect the twisted mass from 
additive mass renormalization. With Eq. (4.2), one can construct the two flavor chiral 
Lagrangian. This is the Gasser-Leutwyler Lagrangian [17] supplemented by chiral and flavor 
symmetry breaking terms proportional to the lattice spacing. The Lagrangian through NLO 
relevant to our work takes the form [124, 153], 

4- = ^tr(a^sa^st) - ^tr(x'ts + stxO - |tr(a^sa^st)2 _ |tr(a^sa,st)tr(a^sa,st) 



h + k 



tr(x'ts + stxOl ' + |tr(a^sa^st)tr(x'ts + st^') 



16 

+ Tytr(5^S5^S"f)tr(i1"S + SU) - Wtr{x^i: + Sl'x')tr(i"fS + SU) 
2 



w 



tr(i1"S + SU) , (4.4) 
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where the LECs, /1-/4 are the SU{2) Gasser-Leutwyler coefficients and the coefficients 
W ^ W and W' are unphysical LECs arising from the explicit chiral symmetry breaking of 
the twisted mass lattice action. The spurion fields are defined as 



X = 2Bo{m + iTsfi) + 2Wob = m + ir^jl + h 
A = 2Wob = b . 



(4.5) 



As discussed in Ref. [124], the vacuum of the theory as written is not aligned with the flavor 
identity but is given at LO by 



with 



(0|S|0) 



M' 



m + 6 + iT^jl 



exp[iujoTs) , 



(rfi + by + fi"^ . 



(4.6) 



(4.7) 



Therefore, to determine the Feynman rules which leave the interactions of the theory the 
most transparent, one expands the Lagrangian around the physical vacuum. Extending this 
analysis to NLO, one finds the vacuum angle shifts to lo = ojq + e where one can determine 
e either by finding the minimum of the potential, as was done in Ref. [124] or by requiring 
the single pion vertices to vanish. 



32: 



e(a;o) = --^bsinuo 



W + 2W'cosuJo — 



One can expand about the physical vacuum by making the replacement 



(4.8) 



^ — sm ^ph Sm; 
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^m = exp(Jtjr3/2), 



(4.9) 



and 



-'ph 



/2i(j) 



One then finds the Lagrangian is given by 




(4.10) 
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b cos uj( WM' + W'b cos uj 
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^^i^M^ph - Kh)) 
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M'f 



If2 



+ 6sinw iWM' + 2H^'6cosw ) tr(Spft + SJ,^) 



(4.11) 



where Ccont. is the continuum SU{2) chiral Lagrangian to NLO. Of particular interest to 
us are the new two, three and four pion interactions which result from the discretization 
errors in the twisted mass Lagrangian. We find, in agreement with Ref. [124] 



C = Ccont. + A£2(ft + A£3<i + A£ 



■'Z4> 



(4.12) 
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where 
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2</> 



coswi^tr(a^(/.a^0) + ]^^M'{u)tT{(t?) + iAM^( 



tr 



V2 



(4.13) 
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smw ^^tr(r30)tr(a^0a^,/.) + ^ ' tr(r3(/>)tr(0^) 
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AM'(w) 
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[tr(0^)] 



AM^ 



2/ 

"^^]^tr(0^) 



(4.14) 



32M/6 , , ^ 
and the mass corrections are given by 



,<9a"; 



cosw — — — tr(9^(/)9^(^)tr(0^ 
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(4.15) 
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AM'(a;) = cosw^ ( VFM' + cos w VF'6 



AMo(tj) 



sin u 



UW'P 
~1^ 



(4.16) 



From this Lagrangian, one can determine the pion masses, decay constants and wave- 
function corrections. One finds the masses are (using the modified dimensional regulariza- 
tion of Ref. [17]) 
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the decay constants are^ 
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and the wave-function correction is 
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(4.17) 
(4.18) 

(4.19) 
(4.20) 



These expressions will be needed to express the scattering in terms of the lattice-physical 
parameters (by lattice-physical, we mean the renormalized mass and decay constant as 
measured from the correlation functions, and not extrapolated to the continuum or infinite 
volume limit). As we discuss in the next section, these interactions lead to three types of new 
contributions to tttt scattering states: there are discretization corrections to the scattering 
parameters, the scattering lengths, effective ranges, etc., which appear in a mild manner as 
those from the Wilson chiral Lagrangian [29]. There are corrections which can potentially 
significantly modify the chiral behavior, arising from the three-pion interactions, and there 
are new corrections which mix different scattering channels, for example, the I = 2,1^ = 
and the / = scattering states. 



^ There is an exact Ward identity one can exploit to compute the charged pion decay constant and avoid 
issues of the axial current renormahzation discussed for example in Ref. [123]. 
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(a) (b) 

Figure 4.1: New unphysical graphs from twisted mass interactions in the t(u)-channel (a) 
and s-channel (6). Fig. (6) can only contribute to I3 = scattering. 

4.3 TTTT Scattering in Twisted Mass xPT 

In this section we calculate corrections to the two-pion scattering channels. We begin 
with the maximally stretched 1 = 2 states, which have the simplest corrections. 

4.3.1 / = 2, /a = ±2 Channels 

There are two types of discretization corrections which modify the / = 2, I3 = ±2 scat- 
tering, those which are similar to the corrections for the Wilson lattice action [29, 155] and 
those which arise from the three-pion interactions, Eq. (4.14) and give rise to new Feyn- 
man diagrams. We will express the scattering parameters in terms of the lattice-physical 
pion mass and decay constant. As was shown in detail, this has dramatic consequences 
on the formula for the scattering parameters in both partially quenched and mixed action 
xPT [79,80,82], such that the extrapolation formulae were free of unphysical counterterms 
through NLO. There is a second benefit to expressing the scattering parameters in lattice- 
physical parameters. This allows one to perform a chiral extrapolation in terms of the ratio 
TTT,^//^, and thus avoid the need for scale setting. This was crucial in allowing the NPLQCD 
Collaboration to make a precision prediction of the 1 = 2 scattering length [110, 111]. 

The simple corrections to the scattering amplitude are determined from AC^^j,, 
Eq. (4.15). The three-pion interactions from Eq. (4.14) lead to new graphs in the scattering 
amplitude, which we depict in Fig. 4.1. The I = 2,1^ = ±2 scattering channels receive 
corrections from Fig. 4.1(a) and its n-channel counterpart. The internal propagating pion 
is a vr'', which, for present lattice actions, is known to be lighter than the charged pions in 
dynamical twisted mass lattice calculations with degenerate light quark masses [119,149]. 

Putting all the corrections together, one finds the scattering amplitude, which we express 
in Minkowski-space, is given by 

T''^' = 7;Lf + Ar2'±2^ (4.21) 

where the discretization corrections are 

Ar-V) = 4^ + c„,.^ s + ^^ (-^ + -^— ) . ,4.22) 

The first two terms arise from Eq. (4.15) as well as from the conversion to the lattice-physical 
parameters. The second two terms arise from Fig. 4.1(a). These terms are formally NNLO. 
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However, depending upon the precision with which the twist angle is tuned, these terms 
may become large and require promotion to lower order. Expanding the NLO contribution 
to the twist angle, Eq. (4.8), one finds 



e^(uj)mt 



'h2 



GAW'b 



sin a; cos u + Oim^^ 



(4.23) 



We can then determine the corrections to the I = 2,1^ = ±2 scattering lengths, for which 
we find 



ArriT^a: 



7=2, ±2/ 



AM'{uj) 



=.2,., a4 



cos(a;) 



^Wbmi {?,2W'Y sin^ uj cos^ w b 



(4.24) 



The first observation we make is that at maximal twist, oj = 7r/2, these leading discretization 
errors exactly cancel through NLO (this is true of the corrections to the scattering amplitude 
and not just the scattering length)^ 



Am^ai-''^\7T/2) = 0. 



(4.25) 



This is independent of the use of lattice-physical parameters, and holds also for the scatter- 
ing length expressed in bare parameters, or any combination of bare and physical. At zero 
twist, cj = 0, our expressions reduce to those of Ref. [29]. Converting /tt — 5- /, our answer 
agrees with that in Ref. [155]. The scattering length at maximal twist is simply given by 
the continuum formula 



rriTra. 



1=2 



-2tt 



m^ 



4vrA 



1 + 



W-T, 



4vr/^ 



3 In 



ml 



/i- 



1-C=1/") 



where the combination of Gasser-Leutwyler coefficients is [89,90] 



7=2 



A{Anf{All + All + ll-ll). 



(4.26) 



(4.27) 



Furthermore, the discretization errors only enter at tree level at this order (when the ex- 
pression is expressed in lattice-physical parameters), and thus at arbitrary twist, the expo- 
nentially suppressed finite volume corrections to Liischer's method are also given by those 
determined in continuum finite volume xPT [113]. 

Returning to the new graphs arising from the three-pion interactions, we can estimate 
the size of the corrections to the scattering amplitude using the known mass splitting be- 
tween the charged and neutral pions [119, 149]. Estimating the splitting with the leading 
correction, Eq. (4.18), and solving for W' from Eq. (4.16), we can estimate the corrections 
to the / = 2, /3 = ±2 scattering length near maximal twist. As a ratio to the LO prediction 
for the scattering length, one finds 



A 7=2,±2 



(w) cot'^uj{AMl^{uj)/ml] 



ml/Sirf^ 



l + AM^(L^)/m2 



(4.28) 



^We have assumed that a suitable definition of the maximal twist angle has been used in the numerical 
lattice computations such that one is not restricted to the regime m, >> b^J^qcDj but rather one is allowed 
m-q > bA%cD [124, 154, 156, 157]. 
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At the lightest mass point calculated in Refs. [119, 149], which corresponds to r7T,,r — 
300 MeV, the pion mass splitting is 

AMniio ~ 7r/2) 

^^ '-^ ~ -0.33 , (4.29 

and therefore one must have cot oj > 0.3 for this term to make more than a 1% correction. 
Therefore, for current twisted mass lattice calculations, corrections to the I = 2,1^ = 
±2 scattering length (and other parameters) should be negligible provided higher order 
corrections are as small as expected. 

4.3.2 I3 = scattering channels 

There are several features which make scattering in the /s = channels more compli- 
cated than in the / = 2, I3 = ±2 channels, most of which stem from the fact that the 
twisted mass lattice action explicitly breaks the full SU{2) symmetry down to U{1), the 
conserved I3 symmetry. The first technical complication is not specific to twisted mass 
calculations, but is simply the need to compute quark disconnected diagrams. The second 
complication stems from the mass splitting of the charged and neutral pions. Generally, one 
determines the scattering phase shift for two particles with the Liischer method [84-86,115], 
by determining the interaction energy 



AE^^ = 2v/p2 + ml - 2m^ . (4.30) 

In the isospin limit, the |2,0) and |0, 0) states (in the [/, /s) basis) are given by 

|2,0) = ^ (Ivr+vr-) + |^-7r+) - 2|7rV)) , 
V6 

|0, 0) = ^ (Ivr+vr-) + |^-vr+) + |^V)) . (4.31) 

v3 

However, given the relatively large mass splitting in current twisted mass lattice calcula- 
tions, Eq. (4.29), the propagating eigenstates will be arbitrarily shifted from the physical 
states, perhaps shifting nearly to the {|vr+vr-), |vr-vr+), [vr'^Tr*')} basis. This would have to 
be disentangled numerically. Even ignoring this issue, which we deem the most serious, and 
working with the continuum {[2,0), |0, 0)} basis, there is a mixing of these states, which 
first appears at NLO as the second operator in Eq. (4.15). Working with the states 

|/.0) = (j^;»>) (4.32) 

one finds 

^^AA. - 9/2 l^__^ 5 J • (4-33) 

Given the correction AMq{uj), one sees this mixing is in fact maximal at maximal twist. 
This is nominally a NLO effect, thus possibly leaving the states mostly unmixed. However, 
a comparison of this term with the LO amplitude of the / = 2, I3 = ±2 scattering, one finds 
close to maximal twist 



-2,0;0,0 



^"TaX; /(32vr) 1 / 1.11 _i.37 



-2vr(m^/4vrA)2 ^2/j2 ^_i.37 1.33 
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(4.34) 
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Figure 4.2: Modified four-point function, (a) consisting of all off-shell graphs. These vertices 
can then be iterated and summed (6), to determine the vrvr interactions. This summation 
gives rise to the Liischer relation, valid below the inelastic threshold. 

For m^//^ = 2, all terms in this scattering matrix are approximately 1/3 the size of the LO 
amplitude. Since we now know that the physical NLO corrections to the 1 = 2 scattering 
length for example, only provide a few percent deviation from the LO term [79,80,110,111], 
we conclude that this NLO operator in fact provides a relatively large contribution to the 
scattering amplitude, and furthermore provides a large mixing term, and thus cannot be 
neglected. This, combined with the problem we mentioned previously, means a coupled 
channel version of Liischer's method of determining the scattering parameters would be 
needed to explore the I3 = scattering channels with twisted mass lattice QCD. 

The A£40 Lagrangian is not the only source of mixing. The three-pion interactions, 
depicted in Fig. 4.1, in the s, t and u channels, will also lead to a mixing of the |2, 0) and 
|0, 0) states, as one can check with an explicit calculation. One may be concerned that the 
new s channel graph will invalidate Liischer's method. This is not the case however, as the 
internal pion propagator is always off-shell, and thus these diagrams do not contribute to 
the power-law volume dependence of the two-particle energy levels. An alternative way to 
understand this is diagrammatically. One can define a modified (momentum dependent) 
four-point function, which is order by order all the diagrams which do not go on-shell below 
the inelastic threshold. We depict this modified vertex in Fig. 4.2(a). These 2PI diagrams 
can then be resummed to all orders to produce the scattering matrix. Fig. 4.2(6). It is this 
resummation that produces the Liischer relation, relating the finite volume scattering to 
the infinite volume scattering parameters [84-86,115]. In this way, one can see that the new 
interactions will not lead to a modification of the structure of the Liischer relation. 

Our last note of caution regards the construction of the interpolating fields. In the 
physical basis, the states which become those of definite isospin in the continuum limit are 
given by Eq. (4.31). However, the interpolating fields are generally constructed with quark 
fields in the twisted basis, with a known definite twist from the physical basis fields. While 
this is also true of the tt'^tt'^ scattering channel, the phase is trivial since there is only one 
term contributing to the |2,2) state. Thus if one were to undertake a calculation of these 
coupled scattering channels, care should be taken in constructing the correct interpolating 
fields. 
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4.4 Conclusions 

In this report, we have detailed vrvr interactions in twisted mass xPT. We have shown 
that through NLO, at maximal twist the corrections to the / = 2,1^ = ±2 scattering 
parameters from discretization errors are identically zero. However, near maximal twist 
there are corrections which can modify the expected chiral behavior which we demonstrated 
by an explicit calculation of the correction to the scattering length. We found however, 
that for the dynamical twisted mass lattice configurations which exist today, the expected 
corrections are negligible. 

The /3 = scattering channels proved to have more significant discretization corrections, 
most notably a mixing term between the |2,0) and |0,0) states which is relatively large. 
In fact, these mixing terms combined with the need for computing quark disconnected dia- 
grams and the expected nonperturbative shift of the twisted mass eigenstates, as discussed 
in Sec. 4.3.2, may make a calculation of these 1^ = scattering channels prohibitively 
complicated. 



47 



Chapter 5 

Meson-Baryon Scattering 
Parameters from Lattice QCD 
with an Isospin Chemical Potential 

5.1 Overview 

As stated in Chapter 4, the determination of hadronic interactions via lattice QCD 
has shown substantial progress in recent years. The main focus on this chapter will involve 
meson-baryon scattering. Several meson-bayon scattering channels are numerically feasible, 
and have been simulated in Ref. [158] . The computation of a different class of hadronic inter- 
actions, however, is hindered by the existence of annihilation diagrams, which, using current 
lattice methods, require all-to-all propagators. While some progress has been made in the 
computation of these propagators [159], the scaling to larger lattices is still prohibitively 
expensive for practical calculations of hadronic interactions. As a result, information on 
pion-nucleon scattering, and KN scattering, a process thought to be important for kaon 
condensation in neutron stars [160-164], has remained elusive to a first principles lattice 
calculation. 

In this work, we propose a novel method for extracting the low-energy constants for scat- 
tering processes with or without annihilation diagrams by employing an isospin chemical 
potential. In terms of lattice calculations, the addition of an isospin chemical potential re- 
sults in a positive definite determinant, which avoids the fermion sign problem that plagues 
other finite density calculations^. An additional feature of an isospin chemical potential 
intrinsic to our analysis is the formation of a pion condensate when the chemical potential 
reaches a critical value [170,171]. This feature has been observed in previous lattice calcula- 
tions [172-177]. In the presence of this dynamically generated pion condensate, propagating 
baryons have their masses shifted by an amount depending on the value of the condensate, 
chemical potential, and a linear combination of LECs, with a dependence that can be cal- 
culated in heavy baryon chiral perturbation theory. As a result, precision measurement of 
baryon masses for several values of the chemical potential coupled with the HBxPT analysis 
gives an avenue for extracting scattering parameters, even for processes with annihilation 
diagrams. 

The organization of this chapter is as follows. In Sec. 5.2, the continuum QCD action 



^For more information and examples, see [165-169]. 
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with an isospin chemical potential is detailed along with the leading-order values for the 
condensates. In Sec. 5.3, we present the two-flavor HBxPT analysis of nucleon and hyperon 
masses in the presence of the pion condensate. Our results show the linear combinations of 
LECs that can be extracted from doing several measurements at different chemical potentials 
above the critical point. In Sec. 5.4, the results of three-flavor HBxPT are derived by 
matching conditions to the two-flavor results. 

5.2 Isospin Chemical Potential 

The primary focus of the next two section will be QCD with two light-quark flavors with 
a finite isospin chemical potential. Using a continuum Minkowski space action, the matter 
part of the QCD Lagrangian is 



c = i) 



1 1 ip+ifiaoY 



'mq + ie-f5Y 



V', 



(5.1) 



where Ip is the SU{3)c gauge covariant derivative, niq = diag(?7i, m) is the quark mass matrix 
in the isospin limit, and ip = (n, d) an isodoublet of quarks. The isospin chemical potential, 
fil, appears as the time-component of a uniform gauge field. The term proportional to e is 
an explicit isospin breaking term. It is included here because the spontaneous breaking of 
isospin is an essential part of our analysis. Because symmetries do not break spontaneously 
at finite volumes, a small explicit breaking is necessary. This explicit breaking should go to 
zero as the size of the lattice goes to infinity.^ 

For vanishing quark masses, m = 0, and vanishing external sources, /U/ = e = 0, the 
QCD action has an SU{2)l <E> SU{2)pi symmetry that is spontaneously broken down to 
SU{2)v- For non- vanishing chemical potential, the QCD action has only a U{1)l iX" U{l)ji 
symmetry associated with I^. This will be broken down to U{l)v by both spontaneous 
chiral symmetry breaking and the explicit symmetry breaking introduced by the quark 
mass. At finite m and finite e, the only continuous symmetry of C is the U{1)b associated 
with baryon number.^ We can write down an effective field theory of QCD that takes into 
account this pattern of spontaneous and explicit symmetry breaking. This effective theory 
is xPT , and at leading order has the form 



£ = — 



< DJJD^'U^ > +2A < M^^U + U'^M > 



(5.3) 



Here the angled brackets denote flavor traces, and the external sources have been included 



^In Euclidean space, the action can be written in the form Se ~ ^{V + mq)'4>, where 

V 






(5.2) 



The Dirac operator D satisfies the condition T^'y^'DT^'ys = 1?^, which ensures that the product of u and d 
fermionic determinants is positive. 

^With both fii and e non-vanishing there is no symmetry that prevents their renormahzation. These 
parameters are only multiplicatively renormalized because, for example, the isospin chemical potential is 
the only term in the action that breaks C ^ t^, and the explicit isospin breaking term is the only term in 
the action that breaks parity. Furthermore because we are interested in the limit of small e, the action will 
be approximately symmetric under the U{l)v associated with t'^. Consequently ni cannot be appreciably 
renormalized. 
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in the terms 

M = s + ip, 
D^U = d^U + i[Y^,U], (5.4) 

where the external vector potential is given by: V^ = fii^S^fi, the pseusdoscalar source is 

2 

given by: p = e^, and the scalar source is just the quark mass, s = iriq. 

The isospin chemical potential favors the presence of up quarks as opposed to down 
quarks. When it is larger than about m.,^, it becomes energetically favorable for the ground 
state to contain positive pions. Assuming a uniform value Uq for the vacuum expectation 
value of U{x), we can minimize the effective potential to determine the vacuum (ground) 
state. Using the standard parametrization of SU{2), namely Uq = exp[ia /^Mn • I^mt], with 
^Mn ■ f^Mn = 1, we find 

TT Jl' \tii\<m^, 

Uo=< r- 21 I I • ^^■^) 

\exp[iaT\, l/U/l > 771^ 

The angle a is determined, at lowest order in the chiral expansion, by the transcendental 

equation 

in \€ 
cos Q = — ^ 7T cot a, (5.6) 

where the pion mass satisfies the Gell-Mann~Oakes-Renner relation, m^ = 2Xm. For small 
values of e, one can determine a using a perturbative expansion of Eq. (5.6). This expansion, 
however, breaks down as one nears the critical value of the chemical potential from above. 
In the condensed phase, the condensates have values 

(ipil;) = f A cos a, 
i{ipT^^^'ilj) = /^Asina. (5.7) 

5.3 Baryons in an isospin chemical potential: SU{2) 

When a hadron propagates in the background of a pion condensate its properties are 
changed due to the interactions of the hadron with the pions in the condensate. Their mass, 
for instance, is shifted by an amount closely related to the forward scattering amplitude of 
pions on the hadron. Unfortunately, the pions in the condensate are not exactly on-shell 
and no model independent relation can be found between pion-hadron scattering amplitudes 
and the mass shift. However, the chiral expansion of these two quantities are related in the 
way described in this section. 

5.3.1 Nucleons 

We can address effects of the isospin chemical potential on nucleons by using the heavy 
nucleon chiral Lagrangian [178]. This effective theory accounts for spontaneous and explicit 
chiral symmetry breaking in the nucleon sector. The theory is written in terms of static 
nucleon fields, and thereby the nucleon mass can be treated in accordance with the chiral 
symmetry breaking scale. 

The leading-order heavy nucleon chiral Lagrangian is written in terms of the nucleon 
doublet, N = {p,n) , and has the form 

£(1) = N^ivD + 2gAS-A)N. (5.8) 
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The vector u^ is the nucleon four- velocity, while the vector, V^, and axial-vector, ^^, fields 
of pions are defined by 



V,, 



A,, 



1 



rt T^L 



ChDjl^L + CRDj^i 



fi^Rj ' 



'- UiDf:^L - e^^d 



(5.9) 



The former appears in the chirally covariant derivative D^, which has the following action 
on the nucleon field 

{D^N)i = d^N^ + {V^)im,. (5.10) 

Appearing in the vector and axial- vector pion fields are ^l, and ^n which arise from the 
definition U = S,lS,r- Under a chiral transformation, {L,R) G SU{2)l <^ SU{2)ji, the pion 
fields transform as: S,l{x) — )■ L^l{x)U{x), and 5,r{x) — )■ U\x)^r{x)R\ where U{x) is the 
matrix entering the transformation of the nucleon field, Ni — )■ U{x)i^Nj. The left-handed 
and right-handed derivatives act according to the rules 



d^Cr 






p- 



(5.11) 



For the case of an external vector field, the left- and right-handed vector sources coincide: 



L„ 



V„. In an isospin chemical potential, the spurious are given the final values 



iL{x) 
Cr{x) 






(5.12) 



so that the vacuum value is Uq = Cq- 

At leading-order using Lagrangian Eq. (5.8), the shift in nucleon mass due to the chem- 
ical potential is the trivial result 



Mn = M 



(0) 

N 



/i/ COSQ- 



(5.13) 



where M^ is the chiral limit value. Beyond leading-order, the nucleon mass receives 
corrections that depend on low-energy constants. These constants are the coefficients of 
terms in the second-order nucleon chiral Lagnrangian. Including relativistic corrections 
with fixed coefficients, the complete'^ second-order Lagrangian is 



£(2) 



ivt 



N 



2M 



[D^,D,] 



WA 



M 



(0) 

N 



\^v A, S ■ D^ + 4ci <M+> 



+4 C2 



9\ 



8M 



(0) 

N . 



(v ■ Af + 4c3^2 + 4c4[S^, S^]A^A^ + Ac^^M. 



N, 



(5.14) 



M^ 



I < M+ >, and M+ = ^X{^lM^]^ + ^rM'i^l)- The last term with 



where A^_ 

coefficient C5 only contributes in the presence of strong isospin breaking. 



*When one considers external sources with non-vanishing field strength tensors, additional terms are 
present. The external sources we consider, however, are uniform in spacetime. 
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Utilizing the second-order nucleon Lagrangian, we arrive at the nucleon mass 

3 / 2 \ 

Mjv = M^' — /i/cosa— + 4ci (mucosa + Aesina) + ( C2 — — ^ + C3 ) /ijsin^ct (5.15) 

The second-order correction is entirely isoscalar, and allows access to the low-energy con- 
stant ci, and the combination C2 + C3. One can vary the quark mass and isospin chemical 
potential to isolate these coefficients from the observed behavior of the nucleon mass in 
lattice QCD.5 

5.3.2 Hyperons 

Recently there has been interest in treating hyperons using SU{2) xPT , effectively 
integrating out the kaon interactions and grouping interactions by strangeness [116,179- 
183]. We consider hyperons in an isospin chemical potential, and determine which low- 
energy constants of SU{2) xPT can be determined from their masses. This is a natural 
approach because the isospin chemical potential produces only pion-hyperon interactions. 

In the stangeness S = 1 sector are the A and S baryons. The A baryon is an isosinglet, 
while the S baryons form an isotriplet. The leading-order Lagrangian in the S = 1 sector 
has the form 

£(1) = aUv ■ dA+ <^^ivD- Aas) S > . (5.16) 

Here we have packaged the S fields into an adjoint matrix 

/I SO s+ \ 
S = h5- 1 .0 h (5-17) 





and use angled brackets to denote flavor traces. The covariant derivative acts as D^T, = 
d^T, + [Vfj_,T,]. The parameter Aas is the SU{2) chiral limit value of the mass splitting, 
Aas = -^s — -^A • In writing the leading-order Lagrangian, we have omitted the axial 
couplings to pions as these terms are not relevant for our computation. Similarly we write 
down only the contributing terms from the second-order S* = 1 Lagrangian 

C = AU [Ac{ <M+> +24 <v-A^> +24 < A^ >] 

+ < S"^S > [4cf <M+> +2c2 <v-A^> +2cf < A^ >] 

-F4c6 <t!^vA><v AT.> +4c7 < SU^ >< A^S > . (5.18) 

Using the first- and second-order Lagrangians, we find the mass of the A baryon 

Ma = Mf^ + a4 {ml cos a + \e sin a) + {4 + 4)f4 sin^ "> (5-19) 

and see that the low-energy constant 4^ ^^'^ the combination 4 + ^3 can be accessed. For 
the S^ baryon, we find the mass 

M^o = M^^ + 4cf {ml cos a + Xe sin a) + (cf + c^)nj sin^ a. (5.20) 



■'One can go further and impose isospin twisted boundary conditions on the quark fields, e.g. xp{x + Lz) = 
exp[i6T ]'ip{x). This has the advantage of introducing non-vanishing spatial components of the vacuum value 
of Afj,, and leads to the ability to isolate more low-energy constants. For the example mentioned, the C3 
coefficient can be determined from the mixing angle between nucleons. 
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At this order in the chiral expansion, the charged S baryons mix. We find the mass eigen- 
values 



M± = Mso + (cf + C7)/i?sin2a±/i/Y'cos2a + (cf + c^)2/x2sin4cK. (5.21) 
The corresponding eigenstates are 



+ 



where the mixing angle can be found from 



COS — \l-i ) - 

2' ' 


- sm — Zj 
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- cos — 2j 

2 ' 



(5.22) 



cosB = =. (5.23) 

' COs2 a + (Cg + Cy )2/i2 gin4 Q, 

Finally because the c^ coefficients should scale with the inverse of the chiral symmetry 
breaking scale, the mixing between S baryons should be small for moderately sized /U/. 

In the strangeness S = 2 sector are the cascades. The case of the cascade is similar to the 
nucleon as the cascades transform as an isospin doublet, H = (5*^,5") . Consequently the 
first- and second-order Lagrangians have the same form as those for the nucleon, Eqs. (5.8) 
and (5.14). One merely replaces the chiral limit mass of the nucleon, Mf^ , with the chiral 

limit mass of the cascade, M4 , the nucleon axial charge, gA with the cascade axial charge, 
ffss- Finally each of the low-energy constants Cj have different values for the cascades; these 
we denote by c^. The cascade masses are thus given by 

M= = Mi — /i/ cos a |-4c'f (m cosa + Aesina) + c? ^%t- + c? Lu/ sin a. (5.24) 

Varying the quark mass and isospin chemical potential allows one to access cf, and the 
combination C2 + C3 . 

5.4 Strange Quark and SU{3) 

We can repeat the above analysis including a third light quark. The isospin chemi- 
cal potential is added proportional to the A^ Gell-Mann matrix, while the explicit isospin 
breaking term we choose in the direction of A^. For a massless strange quark, the QCD 
action with /i/ and e both non-vanishing possesses a U{1)l <^ U{l)ji symmetry associated 
with strangeness. The strange quark mass breaks this symmetry down to U{l)v- This 
symmetry precludes tilting the vacuum value of U in any strange direction. In SU{3), we 
must have 

Co =(=';' ;). (5.25) 

Consequently the non- vanishing values of the condensates are 

^(V^(l-^/3A» = A/^ 

^(V^(2 + ^/3A^)V) = A/2 cos a, 

i{i^X^-f5ip) = A/2 sin a. (5.26) 
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It is now straightforward to determine the effects of the isospin chemical potential on 
the octet baryons. The octet baryons can be packaged into an SU{3) adjoint matrix B, 
given by 

B=\ S- -;^SO + -i=A n . (5.27) 

The relevant term of the leading-order SU{3) Lagrangian is 

C'-^^ =<B^[ivD,B]>, (5.28) 

with Df^B = d^B + [V^, B]. The kinetic term reproduces the isospin-dependent mass shifts 
for the octet baryons found above. 

The contributing terms of the second-order SU{3) chiral Lagrangian are 

£(2) = 2bF < B^[M+,B] > +2hD < B^M+, B} > +2bo < B^B >< M+ > 

+bi < B^[A^, [A^ B]] > +62 < B^[A^, {A^ B}] > +63 < B^A^, {A^ B}} > 
+b4<B^[vA,[vA,B]] >+65 < B^ [v ■ A, {v A, B}] > 
+66 < B^{v -A^iv A, B}} > +67 < B^B >< A^ > +63 < B^B >< v A^ > . 

(5.29) 

Evaluating the masses of the octet baryons, we arrive at results identical to the SU{2) case 
but with SU{3) relations between parameters. The matching of parameters for the baryon 
masses have been given in [180], while those for pion-baryon scattering have been given 
in [182]. Terms proportional to the strange quark mass give rise to shifts of the SU{3) 
chiral limit value of the baryon mass. Such terms are the leading contributions in SU{3) 
to the SU{2) chiral limit values of the baryon masses employed above. As our interest is 
with the /i/ dependence, we do not duplicate the matching relations between the chiral 
limit baryon masses here. Matching of A^-|_ operators in the upper 2x2 block yields the 
relations [180] 



4ci = bD + hF + 2bo, 
M^ = hD + 2bo, 



4 

4cf = 2bD + 2bo, 

4cf = bD-bF + 2bo,. (5.30) 

Matching of operators involving the product of two axial-vector pion fields gives the rela- 
tions [182] 

4(C2+C3) = {bi + b4) + {b2 + b5) + {b3 + be) + 2{bY + b8), 

4(cf + 4) = 1(63 + 66) + 2(67 + 68), 

4(c^+cf) = 4(61 + 64) + 2(67 + 63), 

4(cf+c^) = -2(61 + 64) + 2(63 + 66), 

4(cf + cf) = (61 + 64) -(62 + 65) + (63 + 66) + 2(67 + 63). (5.31) 

From these matching conditions, we see that all three of the LECs for the quark mass 
dependent operators can be determined, along with four combinations of the eight double 
axial-insertion operators. The overdetermination of these parameters, moreover, serves as 
a useful check on SU{3) symmetry. 
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5.5 Discussion 

Above, we show what aspects of meson-baryon scattering can be determined from simu- 
lating lattice QCD with an isospin chemical potential. Our approach allows for the determi- 
nation of the relevant LECs from first principles; and, specifically for pion-nucleon processes, 
avoids the need to evaluate all-to-all propagators. Several attempts at extracting the non- 
derivative pion-nucleon couplings from analysis of experimental data have been attempted in 
the past. In [184], nucleon-nucleon scattering data was used to determine the pion-nucleon 
constants through their effect on two-pion exchange contributions to the nuclear force. Pion- 
nucleon scattering data has also been used for this purpose, for instance [185]. Significant 
uncertainties and discrepancies between different extractions remain. Our method can help 
resolve this situation. The values of the antikaon-nucleon couplings are considerably more 
uncertain, although attempts of extracting them from experimental data have been made, 
see, for example, [186-189]. Our framework can shed light on this issue by testing SU{3) 
symmetry predictions. 

Of particular interest is the coefficient ci, related to the nucleon u-term, which is im- 
portant for addressing dark matter interactions with nuclei, and has been the subject of 
numerous studies. The value estimated in [190] has been challenged [191]. The pion-nucleon 
constants also affect the dependence of the nucleon mass on the quark mass at fourth order 
in the chiral expansion, and this dependence has been used to extract them from lattice 
QCD calculations at different quark masses (for a recent extraction, see [192]). These ex- 
tractions are complicated, however, by the fact that the quark mass dependence of nucleon 
mass is barely compatible with chiral perturbation theory predictions at the current values 
of quark masses [193, 194]. The method proposed here adds another tool that is somewhat 
complementary to the quark mass variation of the nucleon mass. In fact, in order to extract 
the pion-nucleon couplings one needs to know the nucleon mass for several values of the 
quark mass small enough so that chiral perturbation theory is trustable. But, apparently, 
this range might be small [194]. The dependence of the nucleon mass on the isospin chemical 
potential adds another handle that can be used for this purpose. The computational costs 
of changing the quark mass or the chemical potential are comparable. Of course, gauge 
configurations at different quark masses are also useful for other purposes besides nucleon 
mass extraction. The dependence on fij, moreover, can also help in the extraction of other 
observables, making gauge configurations with nj more generally useful, and our proposal 
a realistic option. 
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Chapter 6 

Search for Chiral Fermion Actions 
on Non-Orthogonal Lattices 

6.1 Overview 

One obstacle in lattice QCD is fermion doubling, where the continuum limit describes 
multiple fermions despite only one fermion being attached to each node. While more com- 
plicated discretizations can reduce the amount of doubling, this often comes at the expense 
of exact chiral invariance as shown in the Nielsen-Ninomiya "no-go" theorem [195-197]. In 
the past, attempts were made by Karsten [198] and Wilczek [199] to minimize the doubling 
(only two fermions) allowed by the Nielsen-Ninomiya "no-go" theorem while preserving chi- 
ral symmetry. However, these actions broke additional symmetries leading to the generation 
of non-physical dimension three and four operators in the continuum limit. 

Recently, Creutz [37] and, shortly after, Borigi [200] introduced a non-orthogonal 
graphene-inspired action in four-dimensions in order to achieve both minimal doubling and 
chiral symmetry. Unfortunately, as shown in Ref. [35], this new action suffers from simi- 
lar symmetry breaking issues that plagued the actions proposed by Karsten and Wilczek. 
However, an intriguing realization upon studying the Borici-Creutz action is that if a four- 
dimension action contains the minimum Z5 permutation symmetry, the generation of the 
divergent dimension 3 operators could be prevented in the continuum limit. As a result, 
Ref. [36], explores several non-orthogonal lattices in order to achieve this symmetry, but 
these lattices also lead to other undesireable features. These actions will be discussed 
throughout this note along with the Borigi-Creutz action. 

6.2 Graphene 

Graphene is a single layer of graphite, where the carbon atoms form a two-dimensional 
hexagonal (or honeycomb) lattice. The purpose of this section is to describe some of the most 
basic results of hexagonal lattices that have attracted a great deal of attention throughout 
the condenced matter community. Hexagonal lattices have the intriguing property that 
massless fermions on the sites lead to exact Dirac fermions. 

These lattices are better analyzed by a two sub-lattice description. As shown in Fig. 6.1, 
one sub-lattice (often referred to as the "a-sites" or "1-sites") only communicates with its 
three nearest neighbors which are members of the other sub-lattice (often referred to as 
the "b-sites" or "r-sites"). Another important feature is that these three nearest neighbors 
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Figure 6.1: (left) Graphene lattice composed of two sub-lattices. The "a-site" lattice is 
represented by black circles and the "b-site" lattice is represented by red squares, (right) 
Illustrates the proposal by Creutz [37] to group horizontal neighbors into two-atom sites, 
given by the ellipses. With these sites, a grid can be constructed, which is given by the pink 
lines. 

possess a Z3 permutation symmetry, which means that these lattices are invariant under 
120° rotations. As a result, the addition of the three primary vectors, e", created by the 
three connections to the neareast neighbors yields zero (X^q e" = 0). Without going into 
detail, these features of a 1 -|- 1 dimensional graphene lattice, along with two component 
Dirac spinors, result in chiral Dirac fermions in the massless limit. The question that is of 
interest to this chapter is how to extend this construction to four-dimensions. 

6.3 Graphene-Inspired Lattice Actions 

Several years ago, Creutz in Ref. [37] proposed a viable method of extending this con- 
struction to four dimensional lattice gauge theory. The main idea behind this method was 
to combine the horizontally connected "a" and "b" sites into an individual two-atom "site" 
(as shown in Fig. 6.1 as the ellipses around these two-atom sites). Through this process, the 
connecting segments between these two atom sides form a non-orthogonal grid, which is il- 
lustrated by the pink lines imposed on the lattice in Fig. 6.1. With this formulation, Creutz 
generalized this non-orthogonal grid to four-dimensions while maintaining a construction 
similar to graphene. 

6.4 Borigi-Creutz Action 

The action of this graphene-inspired, non-orthogonal grid in four dimensions proposed 
by Creutz [37], and later Borigi [200] is given by 



Sbc 



lY. 



,1=1 



(6.1) 
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where F^ = (7, ^74) and the four- vectors e" are defined in terms of two parameters B and 
C as 

e' = (1, 1, I, B), e^= ( 1,-1,-1, B), e^ = (-1,-1, 1, B ), 
e^ = (-1, 1,-1, B ), e^ = -( 0, 0, 0,45(7). 

Note, each e° above is a four vector and the a does not refer to the individual components of 
this vector. These e" vectors are generated from this non-orthogonal grid^, which are related 
to the diagonal pink lines in Fig. 6.1. While written in a way that is similar to the naive 
fermion lattice action, there are several key differences in this action. First, the definition of 
F^ contains an ^74, as opposed to the naive fermion action that does not contain this factor 
of i. This leads to an important sign difference in the fourth component of these vectors. 
Second, these e° vectors contain the non-orthogonal structure of this lattice action, which 
lead to non-trivial linear combinations of Dirac matricies. To understand the behavior of 
this lattice action, it is beneficial to write this action in momentum space, which is given 

by 

4 



Sbc = / V'p ^ X] ( sin(p^)e'' • 7 + 574(cos(p^) - C) ) 



V'p- (6.2) 



When written in momentum space, one can immediately confirm that this action has exact 
chiral symmetry. In addition, for B ^ and < C < 1, this action has two poles at 

Pt'^ = P{ 1, 1, 1, 1) # = -j5( 1, 1, 1, 1), 

where cos(p) = C. The existence of only two poles results in the action having minimal 
fermion doubling under the Nielsen-Ninomiya "no-go" theorem [195-197]. 

While the features of exact chiral symmetry and minimal doubling are both convenient 
and attractive features in a lattice action, there are undesirable consequences as well. These 
consequences stem from the action breaking additional discrete symmetries, namely parity 
{4'{'P,Pi) — )■ 74^(— i?,P4)), charge conjugation {'4^{p,Pi) — )■ C-4P"{p,pi)), and time reversal 
{'^{PiPa) ~^ 7574V'(P) ~P4))- These broken symmetries will lead to new unwanted operators 
that are a result of radiative corrections. 

6.4.1 Radiative Corrections 

The presence of additional broken symmetries in lattice actions allow for new operators 
to be generated through radiative corrections. The new unphysical operators due to this 
lattice discretization fall into three categories; relevant, marginal, or irrelevant operators. 
Irrelevant operators, which are operators of mass dimension five or higher, are proportional 
to positive powers of the lattice spacing a and would vanish in the continuum limit, a — )■ 0. 
The other two categories, relevant and marginal, remain in the continuum limit and add 
unphysical terms to the final results that require fine tuning to eliminate. Marginal opera- 
tors, dimension four, can have logarithmic lattice spacing dependence. Relevant operators, 
dimension three or less, are particularly bad since they are proportional to negative powers 
of a and diverge in the continuum limit. 

The exact chirality of the Borigi-Creutz action prevents the relevant operator from the 
Wilson action proportional to a~^VV') but the breaking of parity and time reversal allow to 



^See Ref. [36] for more details. 
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two additional relevant operators, a ^ J2j ^3 C3 , where 

4 4 

O^^^ = 4i5V'74V' = i Yl ^(^^ ■ '>')^' ^?^ = 4i5V^7475V' = ^ Yl ^(^^ " 7)75^'- 

An emphasis should be placed on the fact that fx = 1,2,3,4 and does not include e^ in 
this sum. This illustrates the ^4 permutation symmetry the Borigi-Creutz action pos- 
sesses, which is apparent in this form since these relevant operators are invariant under 
any exchange of the four n values. In addition to these two relevant operators, there are 
additionally ten marginal operators^, which we will not focus on in this note. All these 
operators will be generated unless the action possesses an additional symmetry. 

6.4.2 Additional Symmetry? 

In physical two-dimensional graphene, there exists an additional Z3 rotational symmetry. 
The question that now remains is whether an analogous symmetry exists for the graphene- 
inspired Borigi-Creutz action. Analysis on the subject [36] shows that if an action in four 
dimensions poscesses the minimal Z5 permutation symmetry (or the larger two-operation 
A^ or Ss permutation symmetry), then for a = 1,2,3,4,5, 

5 



E ^" = 0' (6-3) 



a=l 

which results in the relevant operators 

5 5 

0?^ = ij;?(e"-7)V' = 0, of =ij;?(e".7)75V' = 0. 

a=l a=l 

Thus, a lattice action with this minimal symmetry would not have any relevant operators. 
So does the Borici-Creutz action have this minimal Z5 permutation symmetry? To 
explore this question, it is useful to write the action in terms of two two-component spinors, 
one which acts as the "a-site" and one which acts as the "b-site" as illustrated in Fig. 6.1. 
The action written in this way is 



Sbc = ^Yl 



2 

X 



4 

■■Px 



4 
+ Y (Xx-M^ • ""^ ^^ - ^x+„^- e^ Xx) +Xx^-e^(l)x-^x^-e^Xx 

(6.4) 
with S = ((?, —1) and S = {a, 1) and 



V'p = ( ^j , and ijp = {4>p,Xp) ■ 



With this two-component form of the action, examining whether or not this action has this 
symmetry is most easily accomplished by looking at a two-dimensional projection of this 



^See Ref. [35] for more detail on the marginal operators. 
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action and using the graphene picture. In order for any of the terms to display this Z5 
permutation symmetry, the choice of B = l/\/5 and C = 1 is required. The first line of 
Eq. (6.4) is given by the left figure in Fig. 6.2 and shows the desired Z3, nearest neighbor 
behavior in this two dimensional projection (generalizes to the Z5 permutation symmetry in 
four dimensions). However, the next line leads to the right figure in Fig. 6.2, which violates 
Z3 with next-to- nearest neighbor interactions (thus, violates the desired Z5 permutation 
symmetry in four dimensions). Therefore, this action does not have the minimal symmetry 
required to eliminate the relevant operators. 





Figure 6.2: Two-dimensional projection of Borigi-Creutz action with the blue duel basis 
vectors e^ = e'^ — e^. (left) Projection of the first line in of the two-component action shows 
nearest neighbor interactions leading to Z3 rotational symmetry, (right) Projection of the 
second line shows next-to-nearest neighbor interactions breaking Z3. 



6.5 Other Non- Orthogonal Actions 

With the motivation of finding an action with exact chiral symmetry and this minimal 
Z5 permutation symmetry needed to eliminate the relevant operators, Ref. [36] explored 
several non-orthogonal actions which will be reviewed here. 

6.5.1 Modified Borigi-Creutz Action 

With B = l/v5 and C = 1, the first line of Eq. (6.4) obeys the Z5 symmetry and 
the second line breaks this symmetry. A logical extension is to simply delete this last line. 
Therefore, using the same notation as Eq. (6.4), this action with B = l/v5 and C = 1 is 
given by 



"fMBC 



IE 



E 



^,„^S.e^Xx 



/(.X-] 



, S • e^ (/>^) + </.^ S • e^ Xx - Xx S • e^ 



(6.5) 



However, the analysis of this action about the pole p^ ~ yields the behavior ip^i^ ■ k + 
7574^4)^'- This behavior, refered to as mutilated fermions [201,202], is not the desired Dirac 
structure. 
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6.5.2 "Hyperdiamond" Action 

A further modification to the modified Borici-Creutz action to ensure the correct Dirac 
structure in the continuum hmit and preserve at least the minimal Z5 permutation symmetry 
led to the "hyperdiamond" action, which is given by 



E 



4 



11=1 



^X — fJ, ' 



X. 



X+fl 



a 



+ 



: 0- • e Xx 



Xx cr • e^ 



(6.6) 



where 



1 

4' 



\/5, \/5, \/5, 1), e^ 



^(-\/5,-^/5, ^/5, 1), 



-( 0, 



0, 



0, 1). 



\/5,-\/5,-\/5, 1), 



6^ = -(-\/5, \/5,-\/5, 1), 



and CT = (iT, — i), a = (<?, i). The key difference between this action and the modified Borigi- 
Creutz action is the imaginary fourth component of the sigma four-vectors. This alteration 
allows for the action to correctly reproduce the Dirac equation about p^ ~ (which is clear 
in Eq. (6.8)). In addition, the vectors e" satisfy the desired behavior needed for Z5 (or A^ 
or 5*5) permutation symmetry, ^^ e" = and 



e".e/^ 



1, for a = /3 

-1/4, for «//?• 



(6.7) 



While the exact details of the transformation won't be covered here, one can show that this 
action has an A^ permutation symmetry^ invariant under pairs of permutations, which is 
enough symmetry to prevent the relevant operators from being generated. In momentum 
space, the action is given by 



S 



i^p 



4 

i J2 sin(p, 



■M;e^-7 



4 

E 



cos(pp 



e'" + e'^ 



7 75 



^: 



pi 



(6.8) 



with 



i'p = (^] , V'p = {(t)p,Xp) , and 7^ = (_ 







From this form of the action, it is clear that it maintains exact chiral symmetry and cor- 
rectly reproduces the Dirac equation in the continuum limit. Additionally, due to the A^ 
permutation symmetry, relevant operators will not be generated. Unfortunately, this action 
yields excessive fermion doubling, not unlike the naive fermion action. An example of six 



poles in addition to the one at p^ = is pi 
^See Ref. [35] for details. 



-P2 



-P3 =Pa = cos" 



-2/3). 
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6.6 The Wilczek Action 

Karsten [198] and Wilczek [199] pointed out, over twenty years ago, the possibility of 
minimal doubling with an explicit construction. The Wilczek action has the form 



2 

~ 2 ^ [^(^)^'>'4^(^ + •?') + i^{x)h4i^{x - j)] , (6.9) 

with A > 2 as a parameter. The action is diagonal in momentum space 

Sw = f^ ^mD-\k)4^{k), (6.10) 

with the propagator D{k) given by 



Dik) 



y^ ij^ sin A;^ — 174 A V^ (cos kj — 1) 



.11) 



There are exactly two poles of the propagator, namely at the values /c^ = (0, 0, 0, 0), and 

(2) 
kn = (0,0, 0,7r). As this action has been analyzed in detail previously [203], we shall not 

pursue the decomposition of the fields in terms of flavor and chirality. Instead we will focus 

on the discrete symmetries of the action: 

1. Cubic invar iance 

2. CT invar iance 

3. P invar iance 

4. Q invar iance, 

where Q is the time link-reflection symmetry operator. For the Borigi and Creutz actions, 
@ is not a symmetry. 

Given these discrete symmetries, it is straightforward to write down the relevant and 
marginal terms in the Symanzik action. There is one relevant operator 



and several marginal operators 



O3 = tl;ij4il;, 



/n(l) rp Tp 

>-^4 — ^ ^u^ ^u 

C4 = Fii^F^^ 

of =V5D474^ 
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As with anisotropic lattices, the breaking of hypercubic symmetry down to cubic symmetry 
requires a speed of hght tuning. It has been argued in Ref. [203] that the relevant operator 
can be absorbed by a field redefinition, namely 

V'-^expf ^^V, (6.12) 

where C3 is the non-perturbatively determined coefficient of the operator O3 in the Symanzik 

(3) (3) 

action, while c\ that of O^ . The field redefinition, however, modifies the boundary con- 
dition and just moves the fine tuning problem from the bulk of the action to the boundary. 
The physical reality of this term can be ascertained by noticing that it plays the role of an 
imaginary chemical potential. The relevant operator thus cannot be avoided in the case of 
the Wilczek action. 

6.7 Conclusion 

Non-orthogonal lattice actions can be used to enforce desirable features from a lattice 
action. As shown for the Borigi-Creutz action, clever discretizations can enforce exact 
chiral symmetry and minimal doubling. However, upon gaining these benefits, additional 
symmetries are broken, which can (and often will) lead to the generation of relevant or 
marginal operators from radiative corrections. An intricate balance of symmetry is needed 
in a discrete lattice action in order to have exact chiral symmetry, minimal doubling, and 
no relevant operators. At this time, no non-orthogonal lattice action has accomplished this 
task. However, this is by no means a proof that such a lattice action does not exist. The 
possibility still stands that an action with just enough symmetry to rule out these relevant 
operators exists and if found, it would be a very efficient, cheap way to simulate chiral 
fermions. 
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Chapter 7 

Discussion and Outlook 



Effective field theory provides a formalism not only to account for lattice spacing effects 
from lattice QCD, but also can be used extract information indirectly from lattice calcula- 
tions or account for potential issues using new lattice actions. In this work, examples from 
all three of these categories have been presented. 

Lattice QCD is now entering a truly interesting time as a result of increased compu- 
tational power and improved algorithms. In fact, several groups are currently performing 
simulations at the physical pion mass [204,205]. For the first time, physically accurate 
non-perturbative calculations involving nuclear interactions, hadronic spectra, and QCD 
thermodynamics realistically appear to be on the horizon. 

In one sense, extrapolation formula to the physical point from effective field theory are 
not necessary if the calculations are at physical pion masses. While both finite volume and 
lattice spacing effects are just as, if not more, significant as a systematic error, increased 
computing power can reduce these effects substantially. However, from the point of view of 
the effective field theory, there are several interesting questions that can finally be addressed. 
For example, the pion mass is thought of as a small expansion parameter in hadronic effective 
field theories; is this mass small enough to capture all the chiral dynamics? 

An additional question of interest is how to extract parameters from effective field theory 
at the physical point. For example, the sigma parameter in HB^PT is one of great interest 
restricting dark matter models and the KN coupling is important for kaon condensation 
scenarios in neutron stars. One suggestion proposed in this work in Chapter 5, is to employ 
an isospin chemical potential as an additional parameter (or knob) that can be varied to 
help extract these parameters. 

Ultimately, effective field theory and lattice QCD are subjects often intertwined. Effec- 
tive field theory can be used to make clever proposals for lattice QCD measurements and 
alternatively, lattice QCD results can have interesting, unexpected effects that can alter 
the underlying assumptions of an effective field theory. For this reason, the two fields will 
remain intertwined as we enter this new era of hadronic physics. 
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Appendix A 

Scalar Field on the Lattice 



As in Ref. [2], the partition function for a theory with only a scalar field is given by 
(neglecting source terms) 

z= I d\^Y^^'^\ (A.l) 



where d\(\)\ = W^ dcpn represents all possible paths of (p and 

S{<t>) = f d'^xC{<p,d(t>). (A.2) 

Performing the Euclidean rotation, xq — t- ix4, 

e^s{4>) ^ ^-Se{4>) , (A.3) 

The Euclidean correlation function is related to the vacuum expectation value of the time 
ordered fields, 

{<P{xE)<pm = Zf I d[<P]4>{xE)H0)e-''^^'^^ = {0\T[4>{xE)m\0) (A.4) 

The operator (J){xe) can be written as, 

which results in the correlation function, 

{<t>{xE)d^m = (O|,^(O)e-^-^+*P-0(O)|O). (A.6) 

Inserting two complete set of states, 

n 

where Nn represents the proper finite volume normalization. The key result is that the 
calculation of the correlation function is a sum of decaying exponentials whose arguments 
are the energy at each level multiplied by the length in the Euclidean time direction. Since 
the ground state is the smallest energy, if we examine the correlation function at long 
times, the dominant term will be proportional to e"^"^**. Thus, by taking logarithms of 
appropriate ratios at large X4, we can extract the ground state energy, which for a two point 
correlator is the mass. 
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Appendix B 

Three or more flavor xPT at NLO 



The NLO continuum xPT Lagrangian contains operators at 0{p'^), 0{p'^mq), and 
0{mg). The continuum Lagrangian is given by 

Ccont = i.i(tr(a^sta^s))' + L2tr(a^sta,s)tr(a^sta'^s) 
+L3tr(a/,Eta^sa^sta'^s) 

+L4tr(9^S^9^S)tr(m5S"f + Sm^) 
+L5tr(((a^Sta'^S)(m,St + Sm, 

+Le f tr(mgS'l' + Sm^) j + L7 f tr(mgi;'l' - T,mg) j 
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